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Abstract 



Large molecules such as proteins and nucleic acids are crucial for life, yet their primordial origin remains 
a major puzzle. The production of large molecules, as we know it today, requires good catalysts, and the 
only good catalysts we know that can accomplish this task consist of large molecules. Thus the origin 
. of large molecules is a chicken and egg problem in chemistry. Here wc present a mechanism, based on 

, autocatalytic sets (ACSs), that is a possible solution to this problem. We discuss a mathematical model 

^ I describing the population dynamics of molecules in a stylized but prebiotically plausible chemistry. Large 

molecules can be produced in this chemistry by the coalescing of smaller ones, with the smallest molecules, 
the 'food set', being buffered. Some of the reactions can be catalyzed by molecules within the chemistry 
with varying catalytic strengths. Normally the concentrations of large molecules in such a scenario are 
^ ' very small, diminishing exponentially with their size. ACSs, if present in the catalytic network, can focus 

the resources of the system into a sparse set of molecules. ACSs can produce a bistability in the population 
. dynamics and, in particular, steady states wherein the ACS molecules dominate the population. However 



to reach these steady states from initial conditions that contain only the food set typically requires very 
large catalytic strengths, growing exponentially with the size of the catalyst molecule. We present a 
solution to this problem by studying 'nested ACSs', a structure in which a small ACS is connected to 
a larger one and reinforces it. Wc show that when the network contains a cascade of nested ACSs with 
the catalytic strengths of molecules increasing gradually with their size (e.g., as a power law), a sparse 
' subset of molecules including some very large molecules can come to dominate the system. 



Introduction 

One of the puzzles in the origin of life is the question: How did large molecules, which are essential for 
all cells to function, first arise? Macromolecules such as RNA and protein molecules, which contain from 
about a hundred to several thousand monomers, are produced in cells with the help of two crucial catalysts 
(a) the RNA polymerase which reads the genes on DNA molecules and produces the corresponding 
messenger RNA molecules and (b) the ribosome which reads the messenger RNA molecules and produces 
the corresponding protein molecules. These two powerful catalysts, RNA polymerase and ribosome, are 
themselves made up of proteins and RNA molecules, each of which is produced by the process mentioned 
above. When cells produce daughter cells, the latter are already endowed with these catalysts at birth, 
from which they synthesize other molecules. Nowhere in the living world is there a natural process we 
know of that produces macromolecules and that does not itself use macromolecules. Hence the puzzle. 
We expect that the answer to the question lies in the processes that occurred before life originated. 
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The Miller experiment [T] and subsequent work [IHi] were successful in synthesizing monomer building 
blocks of large molecules in simulated prebiotic environments. Those experiments suggested that amino 
acids and nucleotides, monomer building blocks of macromolecules, could be produced on the prebiotic 
earth. Subsequently there has been much experimental work to explore mechanisms that could enhance 
the concentrations of monomers and synthesize long polymers [6l-[9]. While there is interesting progress, 
as yet there is no compelling scenario for the primordial origin of large molecules. 

Meanwhile what has been observed is that catalysis is a fairly ubiquitous property that arises in 
different kinds of molecules and even at small sizes. Organocatalysts [TUHl2j . peptides [T31II1], and 
RNA molecules [T51 - [T7] are known to have catalytic properties. Cofactors play an important role in 
catalyzing metabolic reactions and they (or their evolutionary predecessors) may have had a role in 
prebiotic catalysis [T5] . 

The ubiquity of catalysis motivates the main idea behind the present paper. Here we attempt to 
investigate theoretically, using a mathematical model, whether one can construct a chemical organization 
that produces large molecules from small ones, using the property of catalysis. Apart from the specific 
question of the origin of large molecules the present work is also motivated by a larger question of 
how complex structures and organizations are built incrementally from simpler ones. In systems where 
catalysis is possible an important self-organizing structure that can appear is an autocatalytic set (ACS). 
ACSs were proposed by Eigen (TH], Kauffman [50] and Rossler [3T] and have been used by many authors 
to study various aspects of self-organization, evolution and the origin of metabolism [HH3T], the origin 
of replication [321134] . and the origin and dynamics of protocells [35ti38] . In order to separate the issues, 
the present model only has catalysis and no replication or spatial enclosures; we wish to see what can be 
achieved by catalysis alone. 

Farmer et al [22] , Bagley et al [39] , and Bagley and Farmer [23] proposed and analyzed a model of an 
artificial chemistry in which polymers could form by ligation of shorter polymers through spontaneous 
reactions as well as reactions catalyzed by other polymers in the chemistry. Bagley and Farmer [23] 
analyzed the population dynamics of the molecular species and established some important properties of 
autocatalytic self-organization. When the food set (monomers) were supplied at a fixed input rate and 
the chemistry contained an ACS they showed that in a suitable range of parameters the concentrations 
of the ACS molecules dominated over the rest of the molecules (the background), thereby focusing the 
chemical resources of the system into a small subset of molecules comprising the ACS. However the largest 
polymers in the ACSs they considered had about 15-20 monomers; they did not systematically investigate 
the problems that arise in generating much larger molecules in their chemistry. 

These problems were sharply articulated in the work of Ohtsuki and Nowak [33], in which they 
considered a much simpler model that could be analytically solved. In this model, which they refer to 
as 'symmetric prelife' with a catalyst, they showed that in order for the catalyst to acquire a significant 
concentration in a prebiotic scenario its catalytic strength should be very large, growing exponentially 
with its length. The inference from the model, therefore, was that it is difficult for a large catalyst 
molecule to arise in a prebiotic scenario. 

In this work we consider a model of artificial chemistry similar in structure to that of Bagley and 
Farmer. This model is intermediate in complexity and realism between the model of Bagley and Farmer 
(which is slightly more complex) and model of Ohtsuki and Nowak (which is much simpler). We study 
the dynamics of this model in the presence of ACSs and in particular a structure that we refer to as 
a 'nested ACS' in which a small ACS helps trigger a larger one. We show that this mechanism when 
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iterated across a cascade of nested ACSs avoids the problem of exponentially growing catalyst strengths. 
This mechanism, therefore, provides a possible route to the construction of large molecules in a pre-biotic 
scenario. Apart from these results our work provides an insight, based on the analytic treatment of the 
system under certain approximations as well as numerical work, of certain ACS properties and questions 
such as why ACSs dominate, why nested ACSs work, etc. 

The Model 

The model is specified by describing the set of molecular species, their reactions, and the dynamical rate 
equations for their population dynamics. A special set of molecules, the 'food set', denoted consists of 
small molecular species, / in number, that are presumed to be abundantly present in a prebiotic niche. 
The simplest version of the model (/ = 1) contains only a single monomer species A (or A(l)) whose 
concentration x\ in a well stirred prebiotic region will be assumed to be buffered (constant). The other 
molecules, A(2), A(3),... (dimers, trimers, etc.), whose concentrations are denoted 2:2, X3, . . .. are all made 
through ligation and cleavage reactions of the type A(z) + A(j) ^ K(i + j) with forward (ligation) rate 
constant denoted fc^ and reverse (cleavage) rate constant kf-y The net forward flux of this reaction pair 
is given by vij = {kfjXiXj — k^Xi+j). The rate equations for the system are given by xi = 0, and, for 
n = 2,3,..., 

00 

'^ij ^ ^ '^'^nn ^nXn (1) 

oo 

[k^jXiXj — k.ijXn) + ^ ^ {kin^i+n k^^XiXn) + 2(fc„„a;2n ~ ^nn-'^ri) ~ 4'nXm (2) 

where </>„ represents a loss rate of species n from the region in question. The two terms in the first 
sum represent the formation (respectively, cleavage) of A(n) from (into) smaller molecules. The two 
terms in the second sum and the following bracket represent the cleavage (respectively, formation) of 
larger molecules via reactions that produce (consume) A(n). The stoichiometric factor of 2 before the 
bracket arises because two molecules of A(n) are involved in the corresponding reaction pair. The set 
of parameters kf^ that are non-zero define the set of possible reactions; collectively they define the 
'spontaneous chemistry'('spontaneous' in the sense that the reactions are possible even in the absence of 
catalysts) . A pair of ligation and cleavage reactions can be excluded from the chemistry by setting both 
kfj and A:|^ to zero. The scheme permits chemistries in which some reactions proceed in only one direction 
(ligation or cleavage) by setting only one of kfj and fc^ to zero. However, we will primarily be interested 
in a chemistry in which each reaction is reversible. The existence of the cleavage reactions makes it more 
difficult for the long molecules to survive; thus it is more significant to demonstrate the appearance of 
long molecules in a model in which cleavage reactions are permitted than in one where only the forward 
(ligation) reactions are. 

We consider a simple scheme for catalyzed reactions, assuming that a molecule enhances the rate 
of a reaction that it catalyzes in proportion to its own concentration. Thus, if A(m) is a catalyst of 
the reaction pair A(i) + A(j) ^ A(i + j), then the rate constants of this reaction pair, kf^ and k^^, are 
replaced kfj — > k[j{l + K^^Xm) and kfj — > k(j{l + K^^Xm), where kJ^ is the 'catalytic strength' of the 
catalyst for this reaction pair. The first term in the bracket, unity, represents the spontaneous reaction 
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rate (which is present irrespective of whether the reaction is catalyzed or not), and the second term 
Kl^Xrn represents the enhancement of the reaction rate due to the catalyst. Note that in this scheme 
a catalyst enhances both the forward and reverse reaction rates by the same factor. If a reaction has 
multiple catalysts, n%^Xm is replaced by '^m^m, where the sum runs over all catalysts A(to) of the 
reaction in question. Typically, only a small subset of the spontaneous reactions will be catalyzed. The 
set of catalyzed reactions together with the catalysts and their catalytic strengths will be referred to as 
the 'catalyzed chemistry'. 

When there are / food set (or 'monomer') species a general molecule A is represented as an /-tuple 
of non-negative integers: A = (11,^2, . . . ,i/), where i; is number of monomers of type I contained in A. 
The identity of a molecule in the model is completely determined by the number of monomers of each 
type contained in the molecule; the order in which they appear is irrelevant. Thus the combinatorial 
diversity of distinct compounds containing a total of n monomers (of all types) grows only as a power 
of n (~ nf~^) instead of exponentially (^ /" for strings) if the order had mattered. This simplification 
helps in picturizing the chemistry and significantly reducing the computational power needed to explore 
large values of n. The reaction scheme and rate equations are similar to the 'l-dimensional' version above. 
Details of the general model and explicit examples of rate equations for / — 1 and 2 are discussed in the 
Supporting Appendix SI. 

The main differences between the present model and that of Bagiey and Farmer are (a) a simpler 
representation of molecules (we do not consider molecules as strings), (b) a simpler treatment of catalysis 
(we do not consider intermediate complexes) , and (c) we ignore the effects coming from small populations 
containing a discrete number of molecules. We reproduce the main phenomenon of ACS dominance 
that Bagiey and Farmer observed, but the relative simplicity of the present model allows us to explore 
other phenomena that they do not report about (this includes a multistability in the dynamics and the 
possibility of building large molecules through nested ACSs). 

The main differences with the model of Ohtsuki and Nowak are (a) a much richer spontaneous chem- 
istry of ligation reactions and the inclusion of reverse reactions (which makes an analytical treatment 
more difhcult), and (b) a much more general class of catalyzed chemistries, instead of a single catalyst 
(which allows us to talk of nested ACSs, in particular). With a specific choice of parameters our / = 1 
model reduces exactly to their 'symmetric prelife' model with a catalyst. In spite of greater complexity 
we are able to numerically reproduce their main results in a much more general setting, and also provide 
approximate analytical understanding of the results. 

Autocatalytic Sets (ACSs) 

The dynamics of the above system is particularly interesting when autocatalytic sets (ACSs) are present 
in the catalyzed chemistry. Consider a set S of catalyzed one-way reactions. 'One-way' means that each 
reaction in S is either a ligation or cleavage reaction. Thus the set of reactants and the set of products 
are unambiguously defined for each reaction and the two sets are distinguished. The presence of a given 
ligation or cleavage reaction in S does not mean that its reverse is also necessarily a member of S. Let 
P{S) be the union of sets of products of all reactions in 5, and R{S) the union of sets of reactants of 
all reactions in S. We exclude the food set molecules from both P{S) and R{S). We will refer to the 
set S of catalyzed reactions as an ACS if (a) P{S) includes a catalyst for every reaction in 5, and (b) 
R{S) C P{S). The latter condition implies that all members of R{S) can be produced from the food 
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set by (recursively) applying reactions from within S. An ACS thus ensures the existence of a catalyzed 
pathway, starting from the food set, for the production of each of its products [I9H2T| . Alternative valid 
definitions of an ACS can be given (see [^BT] for one such); the above definition suffices for our present 
purposes and we hope to return to consequences of other kinds of ACSs in the future. Note that if S 
is an ACS, then its extension, S", that additionally includes the reverse of some reactions in S, is also 
trivially an ACS, as in our scheme a catalyst works for both forward and reverse reactions if both exist 
in the chemistry. 

Results 

Spontaneous (uncatalyzed) chemistries 
Nomenclat ure 

Wc first consider the case when none of the reactions is catalyzed, kJ^^ = for all ij pairs. For concreteness 
we first consider spontaneous chemistries that are 'reversible', 'homogeneous' and 'fully connected'. A 
'reversible' chemistry is one for which each allowed reaction is reversible, i.e., kfj ^ <^=> 7^ 0. A 
'homogeneous' chemistry is one in which all the nonzero rate constants are independent of the species 
labels: 0„ = (j) ^or all n, kfj / kfj = kj independent of i and j, and kf^ ^ fc,^ = kr 

independent of i and j. A 'fully connected' chemistry is one in which all possible ligation and cleavage 
reactions are allowed: kf^ ^ and k^^ ^ for all A chemistry is 'connected' if every molecule can 
be produced from the food set in some pathway consisting of a sequence of allowed reactions. In this 
paper we discuss spontaneous chemistries that are reversible and homogeneous. We have checked that 
introducing irreversible reactions and bringing in a small amount of heterogeneity does not change the 
conclusions. Some results for sparse chemistries are discussed later. For homogeneous and fully connected 
chemistries the model has 4 parameters, kf, kr, (j), and the concentration of the monomer, xi = A. 

We explore the model numerically and, to a limited extent, analytically. While the chemistry under 
consideration is infinite, numerical simulations were done by choosing a finite number N for the size of the 
largest molecule in the simulation. In simulating Eq. ([1]) all terms corresponding to reactions in which any 
molecule larger than A(A'^) is produced or consumed were omitted. In principle this introduces another 
parameter, N, an artifact of the simulation. However, one expects that most properties of physical 
interest should become independent of N when N is sufficiently large. Evidence for this is presented 
in Supporting Appendix S2. Our numerical work was mostly done using the CVODE solver library of 
the SUNDIALS (Suite of Nonlinear and Differential/ Algebraic Equation Solvers) package [32], and, for 
smaller N values, using XPPAUT [43] . 

Steady state properties of the spontaneous chemistry: Populations decline exponentially 
with the size of molecules 

Starting from the initial condition in which all concentrations other than the food set are zero (we refer 
to this as the standard initial condition), the concentrations were found to increase monotonically and 
reach a steady state (Fig. [T]A_). Numerically the graph of steady state a;„ versus n on a semi- log plot was 
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found to be approximately a straight line for large n, consistent with the expression 



ce 



(3) 



where, c and A = e^'^ are constants. A, determined by numerically fitting the slope, decreases monoton- 
ically as increases (Fig. [TJ3). 

For (/) = 0, the following exact analytical solution for the steady state concentrations exists for homo- 
geneous and connected uncatalyzed chemistries: 



To see that this is a fixed point, note that when Eq. Q holds, then Vij = kfXiXj — krXi+j = for 
all i,j = 1,2,...; hence the r.h.s. of Eq. ([T]) vanishes (at = 0). Thus A(0 = 0) = kfA/kr- Hence, 
whenever kfA < kr, the steady state concentrations of large molecules are exponentially damped, A < 1. 
When > we do not have an analytic solution. Numerically, we find that A{4> > 0) drops to 



a monotonically decreasing function of kj. and (f>. This corresponds to the intuition that an increased 
ligation rate favours large molecules and an increased cleavage or dissipation rate disfavours them. By 
casting the rate equation in terms of dimensionless variables one can easily see that there are only two 
independent parameters, which may be taken to be k' = kfA/kr and (j)' = (jj/kr whenever kr ^ (for 
details see Supporting Appendix S3). Alternatively when (f> ^ 0, we can take the two dimensionless 
parameters to be and The dependence of A on these two sets of parameters is also shown in 
Supporting Appendix S3. The uncatalyzed chemistry seems to have a global fixed point attractor (all 
initial conditions tested lead to the same steady state). 

Similar results hold when two food sources are present in the system (/ = 2) with buffered concen- 
trations of the monomers (1,0) and (0,1). Simulations are done with all possible reaction and cleavage 
reactions allowed between molecules containing a maximum of TV monomers, all with the same forward 
rate constant kf and reverse rate constant kr and a common dissipation rate (/) for the molecules. A 
steady state concentration profile is shown in Fig. [2l 'Diagonal entries' (ni = 77.2) have higher concentra- 
tions in homogeneous chemistries because there are more reaction pathways to build molecules with equal 
numbers of both monomers than unequal. Since the number of species goes as N'^ /2 and the number 
of reactions as ^ N^, computational limitations require us to work with a smaller N than for / = 1. 
Qualitative conclusions nevertheless appear to be N independent. 

Chemistries with autocatalytic sets 

ACS molecules dominate the population in certain parameter regions 

We now consider chemistries which contain some catalyzed reactions in addition to the spontaneous 
reactions described above. As a specific example to display certain generic properties, we consider the 




(4) 



below 1 even when 



-| — > 1. A is found to be a monotonically increasing function of kf and A, and 
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following catalyzed chemistry: 



A(9) 

A(l) + A(l) ^ A(2) (5a) 
A(5) 

A(2) + A(2) ^ A(4) (5b) 
A(28) 

A(l) + A(4) ^ A(5) (5c) 
A(14) 

A(4) + A(5) ^ A(9) (5d) 
A(37) 

A(5) + A(9) ^ A(14) (5e) 
A(37) 

A(14)+A(14) ^ A(28) (5f) 
A(65) 

A(9)+A(28) ^ A(37) (5g) 
A(14) 

'A(28)+A(37) ^ A(65). (5h) 

Note that this set of reactions constitutes an ACS (which we will refer to as ACS65). If any one 
reaction pair is deleted from the set, it is no longer an ACS. For the moment, for simplicity, we consider 
the case where the catalytic strengths of all the catalyzed reactions are equal ('homogeneous' catalytic 
strengths): Kg'^ = Kg'^ = Kjg"* = Ki'^ = k^7 — '^37'^^ = '^e'd^ — '^14''^'' = all other kJ^ = 0. 

(For clarity, in view of double digit indices, we have introduced a comma between the pair of indices 
in the superscript.) Fig. describes the steady state concentrations, starting from the standard 
initial condition, for the chemistry that contains these eight catalyzed reactions in addition to all the 
reactions of the fully connected spontaneous chemistry. At k = 2.5 x 10^ the ACS product molecules 
dominate over the background (the 'background' being defined as the set of all molecules except the 
ACS product molecules and the food set), in the sense that the ACS molecules have significantly larger 
populations than the background molecules of similar size [23] . There is a fairly sharp threshold value 
of K above which ACS domination appears, as evident from the comparison with the lower curve in Fig. 
[3]A_ drawn for k ~ 2.0 x 10^. Fig. [SjB shows that the steady state background concentrations decline as 
(j) increases, while the ACS concentrations are relatively unaffected in this regime (thus ACS domination 
increases). If catalyzed production pathways from the food set to other molecules are broken somewhere, 
the concentration of the latter molecules declines significantly. This is evident from Fig. [3]3 for which 
only one reaction pair (jSap is deleted from the catalyzed chemistry (which now contains no ACS) while 
others are catalyzed at the same strength as before. 

ACS domination at a sufhciently high catalytic strength also occurs when there is more than one 
monomer. An example with / = 2 is shown in Fig. |4l 
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Understanding why ACS concentrations are large (the k ^ oo limit) 

The above features are generic for a large class of ACSs. It is instructive to consider the k — !> oo limit 
which we discuss analytically. When k is nonzero, the terms in Eq. ([l} corresponding to catalyzed 
reactions get modified. The net flux v of such reaction pairs on the r.h.s. (for brevity we are omitting 
the subscript ij in Vij) is replaced by (1 + Kj2m^"i)v, where the sum over to is a sum over all catalysts 
of the reaction pair. Now let the set S of catalyzed reactions be an ACS. Then, if A(n) G P{S) the r.h.s. 
of in contains at least one such catalyzed term, while if A(n) ^ P{S) in contains no such term. For 
example, for the above mentioned ACS, we have 



X2 


~ k{xqVi_i — 2X5V2.2) + (terms independent of k) 


(6a) 


X4 


— k{x5V2^2 ~ X2SV1A - a;i4W4^5) + (k° terms) 


(6b) 


X5 


~ k{x28Vi,4. - 2:i4W4,5 " xajv^^g) + (k" tcrms) 


(6c) 


Xg 


~ K(a;i4-U4,5 - xgvwg^g - Xe5V<)^28) + («" tcrms) 


(6d) 


ii4 


~ k(x37W5,9 - 2x37Ui4,i4) + (k" terms) 


(6e) 


±28 


~ k{x37Vi4,m - a;65W9,28 " a;i4W28,37) + {^^ terms) 


(6f) 


±37 


~ K(a;65W9,28 - a;i4W28,37) + {k° terms) 


(6g) 


±65 


~ K(a;i4?;28,37) + (k° terms). 


(6h) 



while the rate equations for all other (non ACS) molecules (±3, ig, etc.) have no terms proportional to 
K. In a steady state solution the r.h.s. of Eqs. ([6]) is zero, and to leading order in the k ^ 00 limit 
we must set the coefHcients of k to zero. The coefHcients involve only the ACS fluxes Vij and catalyst 
concentrations. Each coefficient is a sum of terms, and each term is proportional to an ACS flux Vij. 
Thus Vij = for the ACS fluxes provides a steady state solution in the k — > 00 limit. Numerically 
we find that when k is sufficiently high the rate equations converge to this solution starting from the 
standard initial condition. Now Vij = kfXiXj — k^Xi^j, therefore Vij — implies Xi+j = kfXiXj/kr for 
the members of P{S). Since by definition there is a catalyzed pathway from the food set to every ACS 
product, we can recursively express the steady state concentration of every ACS molecule in terms of 
xi = A: Xn = A(fe/A/fcr)"-i. 

It is evident that this argument applies whenever the set S of catalyzed reactions is an ACS; thus for 
every member of P{S), Xn — A{kfA/kr)"~^ is a steady state solution of the rate equations in the limit 
K 00. This is corroborated numerically: in Fig. |3j3 since A ~ kf ~ kr ~ I, all the eight ACS products 
should have x„ = 1 in this limit; the numerical result at k = 3 x 10^ is not too far from this limiting 
analytical value. 

A strong ACS counteracts dissipation 

Recall from Eq. (|4]) and the discussion following it that every molecule in a homogeneous connected 
uncatalyzed chemistry has the steady state concentration x„ = A(kfA/kr)"^~^ when there is no dilution 
flux or dissipation (0 ~ 0), and a smaller concentration when there is dissipation {(j) > 0). We have 
observed above that an ACS with a sufficiently large k can boost the steady state concentrations of its 
members, even when ^ > 0, to the same level. The expression Xn = AikfA/kr)^^^ seems to represent an 
upper limit on the steady state concentration of A(n), which can be approached either when dissipation 
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goes to zero, or. when there is dissipation, by membership of an ACS whose catalytic strength becomes 
very large. 

When the reaction pair A(l) + A(l) ^ A(2) is not catalyzed the production of A(2) takes place at a 
much smaller rate, the spontaneous rate. Therefore its concentration is much smaller, and hence so are 
the concentrations of the larger molecules. 

When A(n) belongs to the background the r.h.s. of i„ contains no term proportional to and all 
the K-independent terms have to be kept, including the 0a;„ term. Thus its steady state concentration 
depends upon (p, and as in the case of the uncatalyzed chemistry, declines more rapidly with n when (/) 
increases. 

Multistability in the ACS dynamics and ACS domination 

The reason for the sudden change in the qualitative character of the steady state profile as k is increased 
is a bistability in the chemical dynamics due to the presence of the ACS. Fig. [5] shows three regions 
in the phase diagram of the system, separated by values and k^^ of k. For < k < (region I), 
the dynamics starting from both the initial conditions mentioned in the figure caption converged to the 
same attractor configuration, which is a fixed point in which the large ACS molecules have a very small 
concentration (the concentration declines exponentially with n). For k^^ < k (region III), again they 
converge to a single attractor, a fixed point in which the ACS molecules have a significant concentration 
which approaches a;„ = A{kfA/kr)"~^ as k — >■ oo. In the range < n < k^^ (region II), they converge 
to two different stable attractors, both fixed points for the ACS under discussion. (We remark that using 
other initial conditions we have found at least one more stable fixed point in a part of region II which 
has intermediate values of xgs, indicating that this system has multistability.) 

This phase structure implies that if we start from the standard initial condition and consider the 
steady state profile to which the system converges for different values of n, we will see a sharp change 
in the steady state profile as k is increased from a value slightly below k^^ to a value slightly above k^^ . 
Below the large ACS molecules will be essentially absent in the steady state, and above k^^ they will 
be present in large numbers and will dominate over the background. 

Therefore, following the nomenclature of Ohtsuki and Nowak [34], who observed a similar bistability 
in their model with a single catalyst, we refer to k^^ as the 'initiation threshold' of the ACS. Similarly 
will be referred to as the 'maintenance threshold' of the ACS, because once the ACS has been initiated, 
K can come down to as low a value as , and the ACS will continue to dominate. 

Bistability in simple ACSs 

In general and k^^ depend upon the other parameters, as well as the topology of the catalyzed and 
spontaneous chemistries. The phase structure is exhibited in more detail for a simpler example in Fig. 
ini where the catalyzed chemistry consists of only two reaction pairs: 



A(4) 



A(l)+A(l) 



A(2) 



(7a) 



A(4) 



A(2)+A(2) 



A(4) 



(7b) 
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which constitute an ACS (called ACS4). This system, investigated numerically using XPPAUT, shows 
bistability. For a fixed <f> the bistability diagram is shown in Fig. |6]4.. The dependence of k^,k^^ on 
(j) is exhibited in Fig. [BB, and on both cj) and kf in Fig. IHP. For a given kf, there is critical value of 
(/>{= (j)) at which the and k,^ curves meet, below which there is no bistability. The locations of the 
phase boundaries, the and curves, depend upon the specific underlying chemistries (catalyzed and 
spontaneous) as well as the ACS topology. The steady state profiles are shown at sample points in the 
phase space in Fig. |Bp. For > it can be seen, that as in the case of the larger ACS discussed earlier, 
if we start from the standard initial condition, the largest molecule of the catalyzed chemistry, here A(4), 
dominates over the background in the steady state only for k > (e-g-, the panel marked 3 in Fig. [SP)- 
In the range < k < k^, it dominates only if we start from initial conditions where it has a large enough 
value to begin with (panel 2b), but not if we start from the standard initial condition (panel 2a). It does 
not dominate for any initial condition \i k < (panel 1). If is below i^, there is a single attractor with 
no significant ACS dominance if n is small (panel 4), or if n is large (panel 5), ACS dominance exists but 
is not very pronounced as the background concentrations are also substantial. 

We remark that while bistability seems to be quite generic in homogeneous chemistries containing 
ACSs, the existence of an ACS does not guarantee that bistability exists somewhere in phase space. 
For example consider the simplest possible chemistry [N = 2) containing only the monomer (which is 
buffered) and the dimer. If we assume that the sole reaction pair A(l) + A(l) ^ A(2) is catalyzed by 
A(2), the catalyzed chemistry is trivially an ACS and the only rate equation is ±2 = kfA^{l + — 
krX2{^ + K.X2) — 41x2- The system can be solved exactly and always goes to a global fixed point attractor 
starting from any initial condition .T2(0) > 0. However, the JV = 3 chemistry defined by the two catalyzed 
reactions 

A(3) 

A(1)+A(1) ^ A(2) (8a) 
A(3) 

A(l)+A(2) ^ A(3), (8b) 

does exhibit bistability at a sufficiently large 4>. Ohtsuki and Nowak [34] had also found a lower limit 
on catalyst size for bistability to exist in their model. Similar results hold for the f ~ 2 case. From our 
simulations a general observation seems to be that bistability is ubiquitous at sufficiently large values 
of (j) hr homogeneous chemistries whenever the smallest catalyst is large enough compared to the food 
set. When it does exist it seems to provide a crisp criterion for 'ACS domination', including 'initiation' 
(k > K^^) and 'maintenance' ((k > k^). 

We must mention that there exists a substantial literature including results on the conditions for 
bistability in chemical systems f44ti46] . It would be interesting to apply some of those results to models 
of the kind being studied here, which involve a large number of molecular species. 

We now make a brief remark about catalyzed chemistries that do not contain an ACS. A particularly 
interesting case arises when the set of catalyzed reactions forms unbroken pathways from the food set 
to higher molecules, but the catalysts are not drawn from the product set of these reactions (i.e., the 
catalyzed chemistry satisfies condition (b) for an ACS but not the condition (a)). In this case we again 
observe the domination of the molecules produced by the catalyzed pathways, but this occurs at even 
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higher catalytic strengths than for ACSs. The reasons are similar to those discussed above in the context 
of the K — !• cx) limit for ACSs (note that in that limit the identity of the catalyst is irrelevant for the 
argument as long as the catalyst has a nonzero concentration). We have also observed multistabilities in 
these chemistries as well as a dependence on similar to the ACS case. In general, catalyzed chemistries 
(with and without ACSs) seem to have a rich and complex phase structure but more comprehensive 
investigations than we have done are needed to nail down the possible range of interesting behaviours. 

A problem for primordial ACSs to produce large molecules: The requirement 
of exponentially large catalytic strength 

A natural initial condition for the origin of life scenario is one where only the food set molecules, and 
perhaps a few other not very large molecules (dimers, trimers, etc.) have nonzero concentrations, while 
the large molecules have zero concentrations. It is from such an initial condition that we would like to see 
the emergence of large molecules through the dynamics. We have seen that in uncatalyzed chemistries, 
the concentrations of the large molecules remain exponentially small (x„ ^ e^'''",7 > 0). In catalyzed 
chemistries, especially in the presence of an ACS, a few specific large molecules produced by the ACS can 
acquire a high population. However, this seems to require a large catalytic strength for the catalysts. For 
example, for ACS65 this happens at k > k^^ = 2226342, starting from the standard initial condition. The 
fact that such a large catalytic strength is needed to produce appreciable concentrations of molecules of 
even moderate length like n = 65 could be a problem for the ACS mechanism to produce large molecules 
in the kind of prebiotic scenario we are considering. In this section we characterize the problem somewhat 
more quantitatively by determining how k^^ depends upon the size n of the catalysts in the ACS. 

As mentioned earlier, the values of , k^^ depend on the topology of the ACS. The topology of the 
ACS includes the set of catalyzed reactions and the assignment of catalysts to each of the catalyzed 
reactions. Define the 'length' L of an ACS as the size of (i.e., the total number of monomers of all types 
in) the largest molecule produced in the ACS. An 'extremal' ACS of length L will be referred to as one 
in which all reactions belonging to the ACS are catalyzed by the same molecule which is the largest 
molecule (of size L) in the ACS. For concreteness, since we are interested in the dependence of k^^ on 
the catalyst size, we consider only extremal ACSs of length L. We assume that the catalyst has the 
same catalytic strength k for all the reactions in the ACS. We wish to determine the bistable region for 
such ACSs and in particular how the values of and depend upon L. These values depend upon 
the precise set of catalyzed reactions constituting the ACS. For illustrative purposes we consider three 
different ways of generating the ACS described under Methods as Algorithm 1, 2 and 3, which generate 
ACSs with different characteristic structure. 

We determine the and k^^ values for ACSs of different values of L numerically. These are plotted 
in Fig. [T] It is evident that increases with L according to a power law ^ L" (with a ranging from 
2.1 to 2.8 for the three algorithms), while k^^ increases exponentially, 

- e"^, (9) 

with p w 0.64 for all the algorithms, a and p depend upon the other parameters. In particular we find 
that p increases with cj), i.e., the catalytic strength needed for large molecules to arise increases faster 
with the size of catalyst at larger values of dissipation. This generalizes, to a much larger class of models. 
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the results of Ohtsuki and Nowak |34| . who found a Uncar dependence of on L and an exponential 
dependence of k^^ . 

The exponential increase of the initiation threshold, k^^ , with L, quantifies the difficulty in using 
ACSs to generate large molecules in the primordial scenario of the type modeled above. This means that 
one needs large molecules with unreasonably high catalytic strengths to exist in the chemistry in order to 
get them to appear with appreciable concentrations starting from physically reasonable primordial initial 
conditions. 



Nested ACSs: Using a small ACS to reinforce a larger one 

We now discuss a mechanism that may overcome the barrier of large catalytic strengths, and may en- 
able large molecules to arise from primordial initial conditions without exponentially increasing catalytic 
strengths. This mechanism relies on the existence of multiple ACSs of different sizes in the catalyzed chem- 
istry, in a topology such that the smaller ACSs reinforce the larger ones, thereby enabling large molecules 
to appear with significant concentrations without exponentially increasing their catalytic strength. 

To illustrate the basic idea we consider the following simple example where the catalyzed chemistry 
contains only two ACSs, one of length three and the other of length eight (which we refer to as ACSS and 
ACSS, respectively), each generated by the Algorithm 2 mentioned above. All reactions of the former are 
catalyzed by A(3) with a catalytic strength K3, and of the latter by A(8) with the catalytic strength Kg. 
Thus the two ACSs are: 

A(3) 

A(l)+A(l) ^ A(2) (10a) 
A(3) 

A(1)-hA(2) ^ A(3), (10b) 

and 

A(8) 

A(l)+A(l) ^ A(2) (11a) 
A(8) 

A(2)+A(2) ^ A(4) (lib) 
A(8) 

A(4)-fA(4) ^ A(8). (lie) 

The catalyzed chemistry consists of the above five catalyzed reaction pairs (we will refer to this catalyzed 
chemistry as ACS3-I-8). The system also exhibits bistability, and the concentration of A(8) in the two 
fixed point attractors is exhibited in Fig. |S]as a function of K3 and Kg- 

When K3 is small the two pictures show the usual bistability of ACS8 along the ks axis. The initiation 
and maintenance thresholds are Kg^ = 1.78 x 10^ and Kg = 1145 given by the location of the boundary 
between the low concentration region (blue, a;g ~ 10"^) and the high concentration region (yellow xg ~ 1) 
along the Kg axis in Figs. [SK and |8j3 respectively. As K3 increases, the initiation threshold of ACSS 
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decreases slowly for a while, then drops sharply near K3 = 141. This value of K3 is the initiation 
threshold of ACS3 when Kg = 0. When K3 exceeds this value, the steady state value of x% is either high 
(yellow, ^ 1) or intermediate (orange, xg ~ 10""^), depending upon the value of Kg. 

The key point is that the initiation threshold of the larger catalyst depends on the catalytic strength 
of the smaller catalyst. The former plummets sharply when the latter approaches the initiation threshold 
of the smaller catalyst, dropping to a much lower value than before (compare the lower limit of the yellow 
region in Fig. [SJV to the left and right of K3 = k|-^ = 141; the value of k|^ plunges several orders of 
magnitude from 1.78 x 10^ at K3 = down to 2178 at K3 ~ 141). Starting from the standard initial 
condition, thus, the larger catalyst can acquire a significant concentration at a much lower value of its 
catalytic strength in the presence of a smaller ACS operating above its initiation threshold than in its 
absence. 

Why a small ACS reinforces a larger one 

We now present an intuitive explanation of the above mentioned property. The argument rests on two 
observations. 

(a) Why the initiation threshold is exponentially large: The first observation attempts to explain 
is so large in the first place. The contribution of a catalyst to the rate of the reaction it catalyzes 

appears through the factor I + ki, where k is the catalytic strength of the catalyst and x its concentration. 
The term unity in the above factor is the relative contribution of the spontaneous (uncatalyzcd) reaction 
rate. If the catalyst is to play a significant role in the reaction, the catalytic contribution to the reaction 
rate should be at least comparable to the spontaneous rate, i.e., kx should be at least comparable to 
unity. As we have seen earlier the concentration of large molecules is typically damped exponentially with 
their size. Therefore the compensating factor k needs to increase exponentially in order for the catalyzed 
reaction rate to be comparable to the spontaneous reaction rate. For concreteness consider the extremal 
ACSs of length L and consider the steady state population of the catalyst A(L) in the low fixed point 
as K is increased. In the spirit of this rough argument one expects that at the initiation threshold the 
term k^^xl should be of order unity. In Fig. IHlwe display this product for different values of L. Though 
there is a secular decreasing trend with L, this product remains of order unity (Fig. [9K) even as the 
individual factors change over several orders of magnitude (Fig. [9j3). This lends numerical support to 
the above explanation for the exponential dependence of k^^ on L. 

(b) Role of the background and spontaneous reactions: The second observation is that when k 
exceeds the initiation threshold for a catalyzed chemistry containing an ACS, not only do the steady 
state concentrations of the ACS product molecules rise by several orders of magnitude, but also those 
of the background molecules rise. As an example compare the two steady state profiles of ACS65 in 
Fig. [SK, which correspond to values of k below and above the initiation threshold. As one goes from 
the lower to the upper curve, the concentration of the ACS members of course increases dramatically 
(as shown by the sharp peaks), but note that the concentrations of other molecules not produced by 
catalyzed reactions also goes up significantly. Thus in the chemistry containing two ACSs (ACS3+8) as 
one moves along the K3 axis in Fig. [SK and crosses the initiation threshold of ACS3 (i.e., K3 exceeds 
Kg^ = 141), the concentration of A(8) (a molecule belonging to the background of ACS3 as its production 
is not catalysed by ACS3) increases from ^ 10~^ (blue region) to ^ (orange region). This increase 
in the concentration of A(8) by a factor of ~ 10"' makes it easier for ACS8 to function and its initiation 
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threshold drops by a corresponding faetor of about 10^ (from ^ 10^ to ^ 10"^). 

This fact highhghts the role of spontaneous reactions in the overall dynamics. The background 
molecules are connected to the ACS through spontaneous reactions, and if it were not for the latter, 
an ACS would not be able to push up the concentrations of its nearby background. We shall refer to a 
structure such as the one described above containing ACSs of different sizes with the smaller ACS feeding 
into the larger one through the spontaneous reactions as a 'nested ACS' structure. 

The role of 'overlapping' catalyzed pathways in nested ACSs 

The above example also serves to highlight some other features of catalyzed chemistries containing mul- 
tiple ACSs. Note that the production pathway of A(8) in ACSS (Eqs. E]) contains one reaction pair in 
common with ACS3, namely the reaction pair A(l) + A(l) ^ A(2). One can consider a situation wherein 
the overlap is greater. E.g., consider the ACS8' defined by 

A(8) 

A(l)+A(l) ^ A(2) (12a) 
A(8) 

A(l)+A(2) ^ A(3) (12b) 
A(8) 

A(2)+A(3) ^ A(5) (12c) 
A(8) 

A(3)+A(5) ^ A(8). (12d) 

Now the set of reactions in ACS3 is a subset of ACS8' (ignoring the catalyst, which is different in the 
two cases). The degree of overlap of the catalyzed reaction sets between a pair of nested ACSs makes 
a difference in the dynamics. Consider, for example, the catalyzed chemistry consisting of ACS3 and 
ACS8', i.e., the set of catalyzed reactions given by Eqs. (jlOp and ([T2|) . which we refer to as ACS3-I-8'. 
Like ACS3-I-8, this chemistry also shows a reduction of Kg^, when K3 exceeds its initiation threshold. 
We find that while at K3 = the value of Kg^ for the two chemistries is not too different (1.6 x 10^ for 
ACS3+8' versus 1.8 x lO'^ for ACS3+8), at Kg 141, reduces to a value 920 in ACS3+8', which is 
less than half of the value 2178 that it reduces to in ACS3+8. Thus a larger degree of overlap between 
the catalyzed reaction sets of nested ACSs causes more effective reinforcement. 

Another example with this behaviour for / = 2 is described in Fig. [TO] In each of the three ACS 
pairs shown in the figure, the smaller ACS, of length 4, is the same, (it will be referred to as ACS(2,2)) 
and is defined by the reactions (each catalyzed by (2,2)) 

(2,2) 

(0,1) + (1,0) ^ (1,1) (13a) 
(2,2) 

(1,1) + (1,1) ^ (2,2). (13b) 
The three larger ACSs, called ACS(5,3)(a), ACS(5,3)(b) and ACS(5,3)(c), respectively, can essentially 
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be determined from the figure. For example, ACS(5,3)(a) consists of the two reaction pairs given by Eqs. 
p3)) . both catalyzed by (5,3) as well as the three reactions 



(5,3) 

(1.0) + (2, 2) ^ (3,2) (14a) 

(5,3) 

(1.1) + (3, 2) ^ (4,3) (14b) 

(5,3) 

(1,0) + (4, 3) ^ (5,3). (14c) 

ACS(5,3)(b) consists of the single reaction pair given by the first of Eqs. (|T3| . catalyzed by (5,3), as well 
as the four reactions 

(5,3) 

(1,0) + (1,1) ^ (2,1) (15a) 
(5,3) 

(1.0) + (2,1) ^ (3,1) (15b) 

(5,3) 

(1.1) + (3,1) ^ (4,2) (15c) 

(5,3) 

(1,1) + (4, 2) ^ (5,3), (15d) 
and ACS(5,3)(c) consists of the five reaction pairs 

(5,3) 

(1,0) + (1,0) ^ (2,0) (16a) 
(5,3) 

(0,1) + (2,0) ^ (2,1) (16b) 
(5,3) 

(1.0) + (2,1) ^ (3,1) (16c) 

(5,3) 

(0,1) + (3,1) ^ (3,2) (16d) 
(5,3) 

(2.1) + (3, 2) ^ (5,3). (16e) 

We consider the population dynamics of chemistries in which the spontaneous part includes all possible 
ligation and cleavage reactions involving molecules with upto iV = 15 monomers with homogeneous rate 
constants kf = kr — 1, (f) ~ and the catalyzed part containing one or more of the above mentioned 
ACSs. When ACS(2,2) is the only ACS present, the system shows bistability with the initiation threshold 
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being 2) = When ACS(5,3)(a), (b) or (c) are the only ACSs present, the initiation thresholds 
for them are 1125197, 1031082, and 1000112, respectively. When ACS(2,2) and one of ACS(5,3) (a), (b) 
or (c) are both present, and the catalytic strength of (2,2) is 552, the initiation thresholds of the three 
larger ACSs reduce to 941, 1256, and 2482, respectively. Again, it is seen that the larger the degree of 
overlap of the two nested ACSs, the more effective is the reinforcement. 

A hierarchy of nested ACSs: A possible route for the appearance of large 
molecules 

The process of nesting discussed above for two ACSs can be extended to multiple levels of ACSs connected 
to each other. Here we discuss sequences of ACSs of increasing size, with the catalyzed reaction set of each 
ACS in the sequence partially or completely contained within the next one, and the catalytic strength of 
molecules increasing with size in a controlled manner. We construct examples of such sequences in which 
large catalyst molecules containing several hundred monomers can acquire significant concentrations 
starting from the standard initial condition, even though all catalysts have moderate catalytic strengths. 

In order to construct a cascade of nested ACSs in which reaction sets of smaller ACSs are completely 
contained in the larger ones (maximal overlap), we used Algorithm 4 described in Methods. This algo- 
rithm produces a cascade of ACSs with g steps (generations), with the fc*'^ generation ACS containing Uk 
new reactions. We studied several catalyzed chemistries containing a cascade of nested ACSs for / = 1 
and 2. One example of each type is presented below; other examples gave qualitatively similar results. 

Dominance of an ACS of length 441 (ACS441) 

For / = 1 we describe a cascade with g = 15 and rii = 1, 712 = ^3 = • ■ • = nis = 2. This catalyzed 
chemistry had 29 product molecules, the largest of which was yl(441) having 441 monomers. The list of 
molecules and reactions is given in Supporting Table S2. The catalytic strength k of each molecule was 
chosen by an explicit length dependent rule 

k(L) = A'L^, (17) 

where K and /? are constants. We describe a simulation with K = 500 and /3 = 1.5. This particular 
rule was chosen to contrast with Eq. ([9]) which characterizes the initiation threshold of an extremal ACS 
of length L. For a value of L such as 441, the exponential function in Eq. ([9]) would have given an 
astronomically large catalytic strength, whereas the much slower growing power law in Eq. (|17p gives 
k(441) = 4.6 X 10^ for the above mentioned values of the constants. Starting from the standard initial 
condition, the steady state concentration profile of this catalyzed chemistry embedded in a fully connected 
spontaneous chemistry with N ~ 800 is shown in Fig. [TTK. This example shows that with the nested 
ACS structure in the catalyzed chemistry, large catalyst molecules can acquire significant concentrations 
starting from an initial condition containing only the monomers, even when catalytic strengths grow quite 
slowly with the length of the catalyst. It is worth noting that product of the catalytic strength of A(441) 
and its steady state concentration (X441 = 0.0077) is about 36000, and this is the factor by which it speeds 
up the reactions it catalyzes over the spontaneous rate. In view of the fact that enzymes containing a 
few hundred amino acids speed up the reactions they catalyze within cells by factors of about 10^ and 
greater, the catalytic efficiency demanded of A(441) does not seem unreasonably high. 
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Eq. (jl7|) . which gives a particular functional form for k{L), is ad- hoc, and, at this stage, merely an 
example given to quantify the level of catalytic strengths that is sufficient for large molecules to arise 
in appreciable concentrations in the prebiotic scenario under consideration if chemistry has the nested 
ACS structure of the kind discussed. One may ask if an even weaker requirement would suffice. We 
have considered smaller values of a (1.2 and 1.0) keeping K fixed, and found that ACS molecules upto 
a particular size (depending on a) do well but that the concentration of larger ACS molecules trails off 
and merges with the background. The size range of ACS molecules that do well can be increased by 
increasing the coefficient K. Since the results depend upon several factors, including the topology of the 
ACS and the uncatalyzed chemistry, a detailed investigation has not been carried out. 

Effect of a 'sparse' spontaneous chemistry 

Fig. [TT] also shows another aspect of ACS dynamics - the relationship between ACS domination and 
the sparseness of the uncatalyzed chemistry. A fully connected uncatalyzed chemistry is one in which all 
possible ligation and cleavage reactions are allowed. A chemistry with average degree k is one in which the 
average number of ligation reactions in which a molecule can be produced is k. In a fully connected / = 1 
chemistry, a molecule of size n can be produced in about n/2 ligation reactions (A(l) + A{n — 1) — > A(n), 
A(2) + A(n ~ 2) ^ A(n), etc.); therefore the average degree of a chemistry containing all molecules 
upto size N is about N/4 (=200 for the chemistry in Fig. [TTK). In Figs. [TTB and [TTIC. we pruned the 
uncatalyzed reaction set to only k ligation reactions per molecule {k = 20 and 2, respectively), randomly 
chosen from all possible ligation reactions producing the molecule. (For molecules too small to have 
k ligation reactions, all ligation reactions were retained. For the ACS molecules the ligation reaction 
producing them in the catalyzed chemistry was included as one of the k uncatalyzed reactions.) Note 
that while we refer to only the ligation reactions and not cleavage reactions for the purpose of defining 
the degree of a molecule, in our simulations all reactions are treated as reverse reactions. That is, for 
every ligation reaction included in the chemistry the reverse (cleavage) reaction is also present. It is 
seen in Fig. [TT] that the increase of sparseness causes the background concentrations to decline. This is 
because there are fewer pathways to produce the background molecules. There is also a larger variation in 
their concentrations because their production reactions are chosen randomly, and background molecules 
produced in reactions that happen to involve the ACS molecules as reactants do better than others. The 
ACS molecules are seen to dominate more strongly over the background in sparser chemistries; this is 
because there are fewer production pathways to the background that drain their concentrations. 

Cascading nested ACSs with / = 2 

An example of a nested ACS with two food sources, ACS(36,28), is presented in Fig. [12] This is also 
generated by Algorithm 4 and has 7 generations with rik = 3 for each generation, the largest molecule 
being (36,28) (the full list of molecules and reactions is given is Supporting Table S3). Again starting 
from the standard initial condition the larger ACS molecules acquire appreciable concentrations with a 
moderate demand on their catalytic strengths. 

As a final example we present in Fig. [13] a cascade of ACSs, named ACS(18,27) after its largest 
molecule, in which smaller ACSs have only a partial overlap with longer ones. This is generated using 
Algorithm 5 (see Methods) and consists of a series of 10 ACSs of increasing length. The detailed list 
of molecules, reactions and catalytic strengths is given in Supporting Table S4. Each successive ACS 
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in the cascade has only a few reactions in common with other ACSs. Unhke in the examples discussed 
above, generated by Algorithm 4, in the present case molecules (except the small molecules) produced 
in the catalyzed chemistry have typically only one or two catalyzed ligation reactions producing them. 
The chemistry also contains a number of catalyzed 'side reactions', which produce molecules that are 
neither catalysts nor reactants in any pathway leading to the largest molecule. In fact there is a 'side 
pathway' consisting of several reactions that may be viewed as 'draining the resources' of the main ACS. 
ACS dominance at moderate catalytic strengths occurs for this chemistry also. The largest k is 50000 for 
the molecule (18,27), and at a steady state population of 0.26 enhances the rate of a reaction by a factor 
of 13000 over the spontaneous rate. This shows that the mechanism outlined by us is not restricted to 
maximally overlapping nested ACSs but is more generic. 

Discussion 

Our work discusses a possible mechanism by which large molecules can arise in a prebiotic scenario. In 
the context of the present model the appearance of large molecules is a natural dynamical consequence of 
chemistry possessing the structure described above - a cascading nested ACS structure (with a not very 
demanding set of catalytic strengths) embedded in a spontaneous chemistry - together with the buffered 
presence of the food set molecules in a well stirred region of space. The mechanism is an incremental 
one: at each step successive step the system is able to access new states made available by the previous 
step while making only an incremental demand on molecular catalytic capabilities. 

The kind of mathematical model we have studied, inspired by the work of Bagley and Farmer, is quite 
abstract; its virtue is the economy of assumptions that go into its structure. The main ingredients are that 
objects can combine with each other in processes or 'reactions' to form other objects and certain objects 
can facilitate certain processes, i.e., 'catalyze' certain reactions. The population dynamics implements 
a simple scheme for how the abundances of the objects would change with time assuming that the 
probability of objects combining is proportional to their abundances. Such a generic scheme while it 
applies in detail to no particular situation allows us to imagine mechanisms at a conceptual level. It 
is significant that in this scheme an ACS can direct the flows towards itself and cause a certain sparse 
subset of objects, including some specific large composite ones, to capture a large fraction of the chemical 
resources. 

At this level of abstraction the model (or a variant with qualitatively similar features) could apply 
to the peptide chemistry as well as an RNA chemistry and to a prebiotic metabolism, as already noted 
by Bagley and Farmer. Indeed it would be equally applicable if a prebiotic environment actually had a 
mixture of ingredients from all these classes of chemistries, a possibility that has been advocated in, e.g., 
[47l|48] . Copley, Smith and Morowitz [47] have proposed a scenario which seeks to explain how the RNA 
world might have originated through a series of incremental steps starting from a primitive metabolism. 
The food sources for this supposed metabolism are CO2, H2, H2S, NII3, etc., in a hydrothermal vent. 
Their scenario envisages multiple stages of increasing complexity which they refer to as (i) the monomer 
stage, in which metabolism, possibly powered by an autocatalytic set such as the reverse TCA cycle, 
produces nucleotides and simple amino acids, (ii) the multimer stage, which produces dimers and small 
cofactors, (iii) the micro-RNA stage, producing of oligonucleotides of length 3-10, (iv) the mini-RNA 
stage, with ll-40mers, followed by (v) the macro- RNA stage, or the RNA world. In their scenario each 
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successive stage produced better catalysts that collectively catalyzed not only their own production from 
the molecules of the previous stage, but also the reactions of the lower stages. This structure is very 
similar to the cascade of nested ACSs that we have discussed. A suitably modified version of our model 
could be constructed to explore the dynamics of this scenario in more detail. At a general level, in the fact 
that the dynamics of our model results in the stable domination, or concentrating, of the large catalysts, 
our mathematical work perhaps lends support to the workability of such a scenario. 

There is another level at which the present model (or its variants) might talk to prebiotic chemistry. 
Morowitz [32] has suggested that the metabolic network itself has a shell like structure which can convert 
simple molecules like CO2, H2, NH3, etc., through "a hierarchy of nested reaction networks involving 
increasing complexity" into purines, pyrimidines, complex cofactors, etc. Reaction sets created by our 
Algorithms 4 and 5 are reminiscent of this picture. Missing from Morowitz's picture is a catalyst assign- 
ment for each reaction from among the molecules in the various shells or from among other catalysts 
accessible prebiotically (e.g., surfaces in hydrothermal vents). It might be worthwhile to attempt to add 
in that information for a more complete scenario and for potential modeling. 

Caveats and future directions 

(1) We have studied the properties of autocatalytic systems, by choosing specific examples of ACS topol- 
ogy and special algorithms for constructing them. This has allowed us to systematically investigate the 
dynamical properties that ACSs offer. We believe that similar dynamical properties would hold for more 
general topological structures than we have considered. Nevertheless the question arises as to whether all 
these structures are very special structures and whether or not they are likely to arise within real chem- 
istry and 'generic' artificial chemistries. ACSs have been shown to be quite generic in a large class of 
randomly constructed artificial chemistries [5(I1[5T| , and a similar analysis could be extended to cascading 
nested ACSs. This would help parametrize or characterize chemistries that would contain such structures 
and those that would not. In this context it would useful to go beyond the simple case we have considered 
in which a molecule is defined by the number of monomers of each type and not their sequence. It may 
also be interesting to look for structures similar to nested ACSs in real metabolic networks using methods 
similar to those in |52| . 

(2) An important related question is one of side reactions (discussed in [531[S3]) which might destroy 
the efficacy of ACSs. In the real chemistry one expects that even if cascading nested ACSs exist, there 
would also exist other catalyzed reactions channelizing the ACS products into pathways leading in other 
directions. Whether substantial populations of large molecules in the nested ACSs can be maintained 
in the presence of such diversion is a question that remains to be systematically investigated. Our last 
example of cascading nested ACSs in fact has several side reactions and it may be noted that large ACSs 
still dominated in that case. We remark that while side reactions can drain resources from ACSs, they 
also help the system to explore new directions in chemical space in an evolutionary scenario. 

(3) We have considered deterministic dynamics in this paper. Stochastic fluctuations are important 
when molecular populations arc small. For a chemistry containing multiple ACSs, Bagley, Farmer and 
Fontana [21] used stochasticity to produce examples of trajectories that differed from each other in the 
sequence of ACSs that came to dominate the reactor. It would be interesting to explore such effects in 
the context of the present model. 

(4) Another simplification we have made is that of homogeneous chemistries, wherein the rate con- 
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stants of all spontaneous reactions have been taken to be the same, and even catalytic strengths, where 
variable, have been taken to be smooth functions of the length. We have checked that introducing a small 
amount of heterogeneity or randomness in the rate constants does not change the qualitative behaviour 
significantly. However, the effect of cranking up the heterogeneity has not been studied. From studies 
of disordered systems in statistical mechanics and condensed matter physics it has become clear that 
such heterogeneity can lead to rugged landscapes, multiple attractors and timescales, and paths that are 
difficult to locate [55) . The dynamics of such systems when they are driven by a non-equilibrium flux 
or buffering of food set molecules is an open question. It is possible that the constraints placed by the 
ruggedness of the landscape will reduce the number of accessible paths. The populating of molecules at 
different levels in a nested hierarchy of ACSs is likely to happen in fits and starts on multiple timescales 
when heterogeneity is accounted for. It is perhaps in such a scenario that one should look for answers to 
the questions raised under (1), (2) and (3) above. 

Methods 

A. Generating extremal ACSs of length L. 

In the following we describe three different algorithms used for generating a set of reactions that provide 
a pathway to produce a molecule of a given length L from the food set. 

Algorithm 1: Incremental, smallest-stepsize, deterministic construction: Each molecule of 
size n [n = 2, 3, . . . , L) is produced in the reaction A(l) + A(n — 1) — > A(n). All such reactions for 
?i = 2, 3, . . . , L are included in the ACS. 

Algorithm 2: Shortest path, top-down, deterministic construction: Start with A(_L). If L is 
even, it is produced in the reaction A{L/2) + A{L/2) A(L). If L is odd. say L = 2m + 1, then A(L) 
is produced in the reaction A{m) + A(m + I) A(L). The same algorithm is used to select a reaction 
for the production of the precursor(s) of A(L) (namely, for A(L/2) if L is even, and for each of A(to) 
and A(m + 1) if L is odd (= 2m + I)), and is iterated for their precursors, etc., until the only reactant 
appearing in the reactions is the food set molecule, A(l). All the production reactions mentioned above 
are included in the ACS. 

Algorithm 3: Incremental, random construction: In this method, starting from the food set 
the set of reactants Rt is sequentially enlarged step by step ( J" = i?o C -Ri C i?2 C . . .) to include larger 
molecules and a reaction chosen randomly until a product of length L is obtained. Rk denotes the set of 
reactants at step k in the algorithm, and Lk the size of the largest molecule in Rk- At step k, pick a pair 
of molecules (say X and Y) randomly from Rk, and determine the product Z formed if they were to be 
ligated {X being the same as Y is allowed). If the size of this product, L{Z), is < or > L, discard 
the pair and choose another pair. If Lk < L{Z) < L, add the ligation reaction X + Y ^ Z to the ACS, 
define Rk+i ~ Rk U {Z}, and iterate the procedure. If L{Z) — L, add the ligation reaction X + Y ^ Z 
to the ACS, and stop. Initially (fc = 0) the reactant set is just the food set, Rq = T . 

To complete the extremal ACS, we assign A(L) as the catalyst of all the reactions generated using 
any of the above algorithms. For the simulations reported in the paper, the reverse of each reaction 
is also included as a reaction catalyzed by the same catalyst. For concreteness we have described the 
Algorithms 1 and 2 for the single monomer case; their generalizations to / = 2 have also been considered 
by us. Algorithm 3, as described above, can be used for any /. 
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B. Generating a cascade of nested ACSs 

Algorithm 4: Incremental, random construction of a sequence of reaction sets with maximal 
overlap: We construct successive sets, or 'generations', g in number, Pi, P2, . . . Pg, of product molecules, 
starting from the food set J-{= Pa). Each generation Pk [k — 1,2,... 5) has a pre-specifed number 
of molecules, nj, (rii, 712, ...Ug need to be specified before running the algorithm). At step k of the 
cascade an ACS Sk is constructed from the previous ACS S^-i by adding reactions between molecules 
belonging to a reactant set Rk consisting of all the products of the previous generations and the food set, 
Rk = U Pi U P2 • • ■ Pk-i- Let Lk denote the size of the largest molecule in Rk- At the beginning of the 
fc"^ step, Pk is empty and the set of reactions in Sk is the same as in Sk-i, except that the catalysts of the 
reactions in 5*^ are not yet assigned. {Sq is the empty set.) To construct P^ and Sk, pick a molecule X at 
random from the previous generation of products Pk-i and another molecule Y at random from Rk, and 
determine the product Z formed if they were to be ligated {X being the same as Y is allowed). If the size 
of this product, L{Z), is < L^, discard the pair and choose another pair. If L{Z) > L^, add the molecule 
Z to Pk and the ligation reaction X + Y ^ Z to Sk- Repeat this procedure until Uk molecules are added 
to Pk and Uk reactions are added to Sk- Assign catalysts to each reaction in Sk (which includes the 
reactions inherited from Sk-i and the new n^. reactions) randomly from Pk- This completes the fc**^ step. 
To get the full cascade with g generations this process is carried out for k — 1,2, .. .g. By construction 
the set of reactions (ignoring the catalyst) of each Sk is fully contained in that of Sk+i but the catalysts 
are different, being drawn from Pk for Sk and Pk+i for Sk+i- The union of the SkS constitutes the set 
of reactions in the catalyzed chemistry. Note that in this catalyzed chemistry reactions have multiple 
catalysts, the multiplicity declining for reactions producing higher generation molecules. 

The size of the food set constrains how large Uk can be. For / = 1, ni = 1 and Pi = {A(2)} as 
the only product one can make from a reaction in the food set is A(l) + A(l) — > A(2), and 71-2 can only 
have values 1 and 2 as the only products one can produce in the second generation (from reactants in 
R2 — {A(1),A(2)}) are A(3) and A(4), etc. Similarly for / = 2, rii can be only 1,2 or 3, as reactions 
among the two food set molecules (1,0) and (0,1) can only produce three molecules (2,0), (1,1) and 
(0,2). 

Algorithm 5: Incremental, random construction of a sequence of reaction sets with partial 
overlap: In this algorithm we decide on a sequence of lengths, Li, L2, - - - , Lg, and generate an extremal 
ACS, denoted Si, of length Li using Algorithm 3 for each i = 1,2,.. .,5. For each i a different random 
number seed is used to initialize Algorithm 3. The union of the Si''s constitutes the set of reactions in 
the catalyzed chemistry. 
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Supporting Information: Appendix SI 
Details of the general model and explicit ex- 
amples for / = 1 and 2 

1 Rate equations in a general case 



When there are / food set species, T — {Ai, A2, . . . , A/}, a general molecule is represented as A = 
(ai,a2, . . . ,a/), where ai (a non-negative integer) is the number of monomers of type A^ contained in 
A. The 'size' or 'length' n of the molecule is defined as the total number of monomers of all types 
in it: n = X^i^i- The food set molecules themselves are represented by the /-tuples (1, 0, 0, . . . , 0), 



A general reaction pair is given by B + C ^ A, where A — (ai, 02, . . . , a/), B = (61, 62, • • ■ , and 
C = (ci, C2, . . . , c^), with ai = bi + Ci. If xa denotes the concentration of the molecule A, the rate 
equations are given by xa = if A e 7^; otherwise, 



where vxy = k^yXxXy — k^yXz is the net forward flux of the reaction pair X + Y ^ Z with forward 
rate constant fc^y and reverse rate constant k^y, 4) a is the loss rate of A, and Qa represents the set 
of unordered pairs of molecules which can combine together to form A (Qa = {(B,C) : bi + Ci = 
ai, for all i = 1,2,...,/}). 

2 Rate equations for a chemistry with / = 1 

For / = 1, instead of using a '1-tuple', we represent the molecules of size n by the notation A(7i) for 
clarity. Following Eq. (jSl.lj) . rate equations for the system are given by ii = and for 71 = 2, 3, . . ., 



(0,1,0,0,...,0), (0,0,. ..,0,1). 




(Sl.l) 



(B,C)eQA 



B.B^A 



00 




(S1.2) 



i^j ,i-\-j—n 



00 




(S1.3) 
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Explicitly, 



±2 = {k[iXiXi - fcfiX2) - {k[2XlX2 - ky2X3) - 2 {k22xl - k22X4:) ~ {k[2X3X2 " ^32-'J^5) - • • • - 022:2 

X3 = (fcf2XiX-2 - k^2^3) - (fcfgXiX's - fcf3X-4) - {k(;jX2X3 ~ k^^^X^) - 2 (/cfgXaXs - fcfgXe) - (fcf3X4X3 - fcf3X7) 

- ... - (^3X3 

±4 = (fc^a;ia;3 - fc^|a;4) + (fcf2a;2a;2 - ^222^4) - {kf^^xix^ - fc^xs) - {k24X2X4 - fc^aJe) - (fc|ia;3X4 - k^^xj) 
- 2 (fc|^a;4 - kf^xs) - (fc|^X5X4 - fc^xg) - ... - 04.T4 



Truncation of the chemistry for numerical simulations: The model has no upper limit on the 
largest molecule that can be produced in the chemistry but for the purposes of numerical simulation we 
assume that the largest molecule that can be produced is of size N. For the finite chemistry we exclude 
all the ligation reactions (and their reverse) that produce molecules of size larger than N. This results 
in iV — 1 coupled ordinary differential equations in / = 1 case, given by: 



(N-n) 

^ ^ {kijXiXj — k^jXn) ~ ^ ^ (kin-^iXn ~ ^in- 



2 {knnxl - k^„X2n) " 0„a(S1.4) 



for n = 2,3, . . . , N . Note that 00 is replaced by {N — n) in the second summation so that the largest 
molecule produced is of length N. 

These ODEs were integrated using CVODE library of SUNDIALS package^ [H] and (for smaller N) 
XPPAUT |43] . The dependence of steady state concentrations on N is discussed separately in Supporting 
Appendix S2. 



2.1 An example with / = 1, TV = 6 

The fully connected spontaneous chemistry has 9 reaction pairs of forward and reverse reactions given in 



the following table: 


Rl: A(l) + A(l) 




A(2) 




R2: A(l) + A(2) 




A(3) 




R3: A(l) + A(3) 


"13 \ 


A(4) 


"13 



R4: A(2) + A(2) 
R5: A(l) + A(4) 
R6: A(2) + A(3) 





22 








X 










^23 


X 


^ . „ 



A(4) 
A(5) 
A(5) 



R7: A(l) + A(5) 
R8: A(2) + A(4) 
R9: A(3) + A(3) 



X 


fcfj, 






^^24 


X 


X 


I.R 






'^33 


X 





A(6) 
A(6) 
A(6) 



^ https: / / computation. Unl. gov /case/ sundials / main, html 
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From Eq. (|S1.4[) . the 5 equations defining the dynamics are: 



-(1)2X2, (SI. 5a) 

is ^{k^2^iX2 - fcfjXs) - (fc^xixa - fcf3.T4) - (fc|3.T2a;3 - fefa^'s) - '^'{k^zxl - k§^XQ) - 03X3, (SI. 5b) 

±4 =(fc;{^a;ia;3 — fc;f|a;4) + (fc^X2 — ^222^4) ~ {kf^xiXA, — ^14X5) — (fc|^a;2a;4 — fc^xg) — ^4X4, (SI. 5c) 

i.5 =(fc^xiX4 - fc^xs) + (fc|^X2X3 - fc^xs) - (fcfgXix.g - fc^xe) - 0.52:5, (S1.5d) 

X6 =(fcf5XiX5 - fcfgxg) + (^^4X2X4 - ^^4X5) + (fcf3X3 - fcfgxe) - (jy&XQ. (S1.5e) 

2.2 An example of inclusion of catalyst 

Say reaction pairs Rl and R4 of the previous example are catalyzed by A(4) with catalytic strengths 
and respectively. Then the equations for X2 and X4 get modified to 

X2 =(1 + K4^X4)(fcf"j^X^ — fc;Pj^X2) — {k^2^1^2 ~ ^122^3) ~ 2(1 + K|^X4)(fc^X2 — fc^X4) — (fc232:22;3 — ^232^5) 

- (fc|^X2X4 - fc^xg) - 02a;2, (SI. 6a) 

X4 =(fcf3XlX3 - ^0^3X4) + (1 + K4^X4)(fc^2^2 ~ ^22^:4) - (fcf^XiX4 — fc^^X5) — (fc|^X2X4 — fc^Xg) 

— 04X4. (SI. 6b) 

The equations for X3, X5 and xg remain as before. This chemistry, but with N extended to 15 is 
discussed in the context of bistability in Fig. 6 of main text. There all the spontaneous rate constants 
and xi have been chosen to be unity, all the 0„'s are equal and = = k. 



3 Rate equations for a chemistry with / = 2 



For a / = 2 chemistry, J- = {Ai, A2}, a general molecule is given by A = (oi, 02), and the monomers are 
given by 2-tuples (1,0) and (0,1). There are {n + 1) molecules of length n, given by 2-tuples: (n, 0), (71 — 
1, 1), . . . , (1,71. — 1), (0,71,). A finite chemistry, i.e., a chemistry in which the largest molecules that can 
be produced is of length N, has a total of i^+^'>(^+^) _ i molecules (including monomers), and hence 

(Af+l)(Ar+2) 
2 



3 rate equations. 



3.1 An example with / = 2, iV = 3 

The fully connected spontaneous chemistry has 9 reaction pairs given in the following table: 



Rl: (0,1) + (0,1) ^ 
R2: (0,1) + (1,0) ^ 
R3: (1,0) + (1,0) ^ 



(0,2) 

(1,1) 
(2,0) 



R4: (0,1) + (0,2) 
R5: (0,1) + (1,1) 
R6: (1,0) + (0,2) 









X 


X 






X 


X 



(0,3) 
(1,2) 
(1,2) 



R7: (0,1) + (2,0) 
R8: (1,0) + (1,1) 
R9: (1,0) + (2,0) 



X 






X 


X 






X 


X 



(2,1) 
(2,1) 
(3,0) 



Using Eq. (jSl.ip suitably modified for a finite chemistry (as for / = 1 case in Section [2]), we can write 
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the following rate equations: 

i(o,2) =i^Ri - vm - vne - 4>{o,2)^{o,2) (SI. 7a) 

=i'R2 - VR5 - VRs - 0(i,i)2;(i,i) (SI. 7b) 

i(2,0) =1'R3 - - VR9 - 4>{2.a)X{2,0) (SI. 7c) 

i(0,3) =WR4 - 0(0,3)2^(0,3) (S1.7d) 

i{i,2) =WR5 + - (/>(!, 2)a;(i, 2) (S1.7e) 

a;(2,l) =^'R7 + WR8 - '^(2,i)a;(2,i) (S1.7f) 

i(3,0) =t'R9 - 0(3,0)2^(3,0)- (Sl-7g) 

WRi is the net forward flux of the reaction pair Ri (as given in the table above). Eqs. (jS1.7p can be 
expanded as to get the equations in x's. For example, Eq. (jS1.7b[) becomes 

i(i,i) = (fcR22;(o,i)a;(i,o) - fcR2a;(i,i)) - (fc^5a;(i^o)a;(i,i) - fcR5X(2,i)) - {kn8X^o,i)X(is) ~ kisX(h2)) 

- 0(i,i)a;(i,i) (SI. 8) 

where, and are the forward and reverse rate constant of reaction Hi, respectively. 

When the catalysts are included, the reactions that involve catalysts are changed in the same way as 
in / = 1 case (see Section 12.21 of this appendix) . 
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Supporting Information: Appendix S2 
A^- independence of results at large 

While the chemistry considered in this model is an infinite one, numerical exploration requires us to work 
with a finite set of molecules that participate in the chemistry. This introduces an additional parameter 
N, the size of the largest molecule produced in the chemistry, in the model (for details see Appendix 
SI). Here we present evidence showing that certain important properties of the model become essentially 
A^-independent at large enough TV. 

1 Uncatalyzed chemistry 

The steady state concentrations for an uncatalyzed chemistry can be approximated by an exponential 
hmction, Xn = ce^'*'" = cA", where c and A = are constants (see Eq. (3) in the main text). We 
calculate A from the slope of the straight line fit on the plot of the log of the steady state concentrations, 
Xn, versus length, n. We found that for sufficiently large N, A becomes independent of N; see Fig. IS2.1I 
When A < 1, the concentrations of the large molecules are small. It is evident that if the concentrations 
of the large molecules produced in the chemistry are so small that their contributions to the dynamics 
of the (relatively) smaller molecules is negligible, then, including molecules larger than those already 
considered in such a chemistry would cause no difference to the results. 

2 Catalyzed chemistry 

For chemistries with catalyzed reactions, including ACSs, we find that for sufficiently large the nu- 
merical results are independent of N. In Fig. IS2.2I we show how the steady state concentrations in the 
chemistry that includes ACS65 (defined by Eq. (5) in main text) depend upon N. 

The results for the chemistry with / = 2 were also found to be A^- independent (for sufficiently large 
values of A^). We find that for / = 2 the A^- independence is reached at much smaller values of A^ than 
for / = 1. 
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Supporting Information: Appendix S3 
Dimensionless rate equations and dependence 
of A on dimensionless parameters 

Eq. (2) in the main text can be cast in dimensionless form by introducing a concentration scale lo and a 
time scale t. We discuss below the case of a homogeneous spontaneous chemistry. 

Define dimensionless quantities 7i„ = .T„/aj, t' = t/r, k'^ = kfUJT, fc', = krT and (/)' = tpr. Then 

. _ dxn _ d{unUj) _ oj du^ . , 

" ~ d{t'T) ~ T dt" ^^'^ ' 

and Eq. (2) in the main text becomes 

~~Jir = ( —{ujUi){uJUj) - —{UJU„)] ~ ( ^{uJUi){ujUn) - —{uJUi+ri) 

T dt' ^ \UJT T J \UJT T 

_2 f ^(c^u„)' - ^M2„)) - ^{^Un) (S3.2) 



d 



Without loss of generality one can choose uj = xi = A and (whenever fc^ 0) r = l/fc^. Then the 
dimensionless concentration variables satisfy the same equations as before, but with A = k,, = 1. There 
are now only two independent dimensionless parameters, fc^ and </>'. The dependence on all 4 parameters 
can be recovered at the end by replacing u„ by x„/A, t' by tfcr, fc^ by kfA/kr and (j)' by (j)/kr. The 
behaviour of A as a function oi kf and (f) (keeping A = k,- ^ 1) is shown in Fig. IS3.1I 

One may also choose r = 1/0 (whenever 0^0). In that case the two independent dimensionless 
parameters will be fc^ and fc,',. The behaviour of A as a function of kf and kr (keeping A = (j) = 1) is 
shown in Fig. IS3.2I 

When kr = and = both the above choices for scaling fail. In this case the system has no steady 
state for any finite iV as xat is always positive. 
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Supporting Information: Table SI 

List of reactions and their catalysts in ACS (8,10) 

(referred in Fig. 4 of main text) 



Reaction 


Catalyst 


(0,1) + (1,0) ^(1,1) 


(1,3) 


(1,0) + (1,0) ^(2,0) 


(3,6) 


(0,1) + (1,1)^(1,2) 


(5,6) 


(0,1) + (1,2) ^(1,3) 


(8,10) 


(1,1) + (1,2) ^(2,3) 


(1,3) 


(1,3) + (2,3) ^(3,6) 


(1,2) 


(2,0) + (3, 6) ^ (5,6) 


(2,3) 


(2,3) + (5,6) ^(7,9) 


(3,6) 


(1,1) + (7,9)^(8,10) 


(2,3) 



Note that this catalyzed chemistry forms an ACS. Each catalyst in the above list gets produced by a 
reaction of this chemistry. Also, every reactant is either a member of the set T or gets produced in this 
chemistry. 



32 



Supporting Information: Table S2 

List of reactions and their catalysts in ACS441 

(referred in Fig. 11 of main text) 

The tabic lists all the reactions with their respective catalysts in the example of a catalyzed chemistry, 
quoted in the main text, containing a cascade of nested ACSs for f ~ I generated using Algorithm 4. 
The steady state concentrations for this chemistry are displayed in Fig. 11. This chemistry was generated 
with g = 15 and ni = l,n2 = ^ . . . ^ ni5 = 2. 

Note that a general molecule for the / = 1 case, represented by A(n) in the main text, has been 
represented here by n for brevity. 

The molecules in various generations are as follows: Pq = {1}, Pi = {2}, P2 = {3,4}, P3 ~ {6,7}, 
P4 = {10, 12}, P5 = {17, 18}, Pe = {20, 24}, P7 = {25, 48}, Pg = {52, 66}, Pg = {67, 69}, Pw = {77, 84}, 
Pii = {144,168}, P12 = {221,288}, P13 = {305,336}, Pu = {372,389}, P15 = {397,441}. 

The catalyst for a reaction listed under generation P^ is added at step k of algorithm. It is apparent 
from the reaction table that the ACSs are maximally overlapping, i.e., any ACS of generation k contains 
all the reactions of generation k — I. 



Reaction 


Catalyst added in generation 


Pi 


P2 


P3 


P4 


P5 


Pc, 


Pr 


Ps 


Po 


Pio 


Pii 


Pi 2 


Pl3 


Pl4 


Pi 5 


1 + 1^2 


2 


3 


6 


10 


17 


24 


48 


52 


69 


84 


144 


288 


336 


372 


397 


1 + 2^3 




3 


7 


10 


18 


20 


25 


52 


67 


77 


168 


221 


305 


372 


397 


2 + 2^4 




4 


7 


10 


18 


20 


25 


66 


67 


77 


168 


221 


336 


372 


397 


3 + 3^6 






6 


12 


17 


24 


25 


66 


67 


84 


144 


221 


336 


372 


441 


3 + 4^7 






7 


12 


18 


24 


48 


52 


67 


84 


144 


221 


336 


372 


397 


4 + 6 ^ 10 








10 


18 


24 


48 


52 


69 


77 


168 


288 


336 


389 


397 


6 + 6 ^ 12 








12 


17 


20 


48 


66 


69 


77 


144 


221 


305 


372 


441 


7 + 10 ^ 17 










17 


24 


25 


52 


69 


84 


168 


221 


336 


389 


441 


6 + 12 ^ 18 










18 


20 


48 


66 


69 


77 


168 


288 


336 


372 


441 


2 + 18 ^ 20 












20 


48 


52 


69 


84 


144 


221 


336 


389 


441 


7 + 17 ^ 24 












20 


25 


66 


69 


77 


168 


221 


336 


372 


397 


1 + 24 ^ 25 














48 


52 


69 


84 


144 


288 


305 


372 


441 


24 + 24 ^ 48 














48 


66 


67 


77 


144 


221 


305 


389 


441 


4 + 48 ^ 52 
















66 


69 


84 


168 


221 


336 


372 


397 


18 + 48 ^ 66 
















66 


69 


84 


144 


288 


336 


389 


441 


1 + 66 ^ 67 


















67 


77 


168 


221 


336 


389 


441 


3 + 66 ^ 69 


















67 


84 


168 


221 


305 


372 


441 


10 + 67 ^ 77 




















77 


168 


288 


305 


372 


441 


17 + 67 ^ 84 




















77 


144 


221 


305 


389 


397 


67 + 77 ^ 144 






















144 


288 


305 


389 


397 


84 + 84 ^ 168 






















144 


221 


336 


389 


441 


77 + 144 ^ 221 
























221 


305 


372 


397 


144 + 144 ^ 288 
























221 


336 


372 


441 


17 + 288 ^ 305 


























305 


372 


441 


48 + 288 ^ 336 


























336 


389 


397 


67 + 305 ^ 372 




























389 


397 


84 + 305 ^ 389 




























389 


397 


25 + 372 ^ 397 






























441 


69 + 372 ^ 441 






























397 
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Supporting Information: Table S3 

List of reactions and their catalysts in ACS (36,28 

(referred in Fig. 12 of main text) 

The tabic lists all the reactions with their respective catalysts in the example of a catalyzed chemistry, 
quoted in the main text, containing a cascade of nested ACSs for / = 2 generated using Algorithm 4. 
The steady state concentrations for this chemistry are displayed in Fig. 12. This chemistry was generated 
with g = 7 and rifc = 3. 

The molecules in various generations are as follows: Pq = {(1, 0), (0, 1)}, Pi = {(1, 1), (0, 2), (2, 0)}, 
P2 = {(1, 3), (2, 2), (3, 0)}, Ps = {(2, 6), (2, 3), (5, 2)}, P4 = {(4, 6), (7, 4), (7, 8)}, P5 = {(8, 12), (6, 12), (14, 8)}, 
Pe = {(14, 24), (22, 20), (11, 12)}, P7 = {(29, 28), (36, 28), (24, 26)}. 

The catalyst for a reaction listed under generation P^ is added at step k of algorithm. It is apparent 
from the reaction table that the ACSs are maximally overlapping, i.e., any ACS of generation k contains 
all the reactions of generation k ~ 1. 



Reaction 


Catalyst added in generation 


Pi 


P2 


P3 


P4 


Ps 


Pe 


P7 


(0,1) + (0,1) -(0,2) 


(1,1) 


(1,3) 


(2,3) 


(7,4) 


(8, 12) 


(11,12) 


(29,28) 


(0,1) + (1,0) ^(1,1) 


(0,2) 


(1,3) 


(2,6) 


(4,6) 


(14,8) 


(22,20) 


(36,28) 


(1,0) + (1,0) ^(2,0) 


(1,1) 


(3,0) 


(5,2) 


(7,4) 


(14,8) 


(14,24) 


(29,28) 


(1,0) + (2,0) -(3,0) 




(2,2) 


(5,2) 


(4,6) 


(14,8) 


(14,24) 


(24, 26) 


(1,1) + (0,2) ^(1,3) 




(1,3) 


(2,6) 


(7,8) 


(8, 12) 


(14,24) 


(24, 26) 


(0,2) + (2,0) -(2, 2) 




(1,3) 


(2,6) 


(7,8) 


(6, 12) 


(22,20) 


(29,28) 


(0,1) + (2, 2) ^(2, 3) 






(2,6) 


(7,4) 


(14,8) 


(22,20) 


(24, 26) 


(3,0) + (2, 2) ^(5, 2) 






(2,3) 


(7,8) 


(6, 12) 


(11,12) 


(29,28) 


(1,3) + (1,3) -(2,6) 






(5,2) 


(7,8) 


(14,8) 


(22,20) 


(24, 26) 


(2,0) + (2, 6) -(4, 6) 








(4,6) 


(14,8) 


(22,20) 


(36,28) 


(2, 2) + (5, 2) ^(7, 4) 








(7,4) 


(6, 12) 


(14,24) 


(36,28) 


(5, 2) + (2, 6) -(7, 8) 








(7,4) 


(14,8) 


(14,24) 


(36,28) 


(2, 6) + (4, 6) ^(6, 12) 










(14, 8) 


(11,12) 


(29,28) 


(4, 6) + (4, 6) ^(8, 12) 










(14, 8) 


(22, 20) 


(36,28) 


(7, 4) + (7, 4) ^(14, 8) 










(8, 12) 


(11,12) 


(36,28) 


(3,0) + (8, 12) ^(11, 12) 












(22, 20) 


(36,28) 


(8, 12) + (6, 12) ^ (14,24) 












(22,20) 


(24,26) 


(14, 8) + (8, 12) ^ (22,20) 












(11,12) 


(29,28) 


(2,6) + (22,20) — (24,26) 














(29,28) 


(7,8) + (22,20) — (29,28) 














(29,28) 


(14,8) + (22,20) — (36,28) 














(29,28) 
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Supporting Information: Table S4 

List of reactions and their catalysts in ACS (18, 

(referred in Fig. 13 of main text) 

The table lists all the reactions with their respective catalysts and catalytic strengths in the example of 
a catalyzed chemistry, quoted in the main text, containing a cascade of partially overlapping ACSs for 
f ~ 2 generated using Algorithm 5. The steady state concentrations for this chemistry are displayed in 
Fig. 13. 

The catalyzed chemistry contains 10 generations of ACSs of lengths 3, 6, 10, 15, 19, 25, 30, 35, 40, 
and 45. 

Table S4a: List of reactions in the ACSs of different length and their catalysts. ACSs 
of increasing length (using Algorithm 5) are added in the chemistry. The length of the ACS and the 
catalytic strength of the catalyst arc mentioned in the table. 



Reaction 


Catalyst 


Generation 1 of length 3 (k = 1000) 


(1,0) + (1,0) ^ (2,0) 


(3,0) 


(2,0) + (1,0) ^ (3,0) 


(3,0) 


Generation 2 of length 6 (k = 2000) 


(1,0) + (1,0) ^ (2,0) 


(4,2) 


(0,1) + (2,0) ^ (2,1) 


(4,2) 


(2,1) + (0,1) ^ (2,2) 


(4,2) 


(2,1) + (2,1) ^ (4,2) 


(4,2) 


Generation 3 of length 10 {k = 4000) 


(1,0) + (0,1) ^ (1,1) 


(4,6) 


(1,1) + (1,0) ^ (2,1) 


(4,6) 


(1,1) + (1,1) ^(2, 2) 


(4,6) 


(0,1) + (2, 2) ^ (2,3) 


(4,6) 


(2,3) + (2,3) ^ (4,6) 


(4,6) 


Generation 4 of length 15 (k = 7000) 


(0,1) + (0,1) ^ (0,2) 


(7,8) 


(0,2) + (1,0) ^ (1,2) 


(7,8) 


(1,2) + (1,0) ^ (2,2) 


(7,8) 


(2,2) + (2,2) ^ (4,4) 


(7,8) 


(1,2) + (4, 4) ^ (5,6) 


(7,8) 


(5,6) + (1,2) ^ (6,8) 


(7,8) 


(5,6) + (2,2) ^ (7,8) 


(7,8) 


Generation 5 of length 19 (k = 10000) 


(0,1) + (1,0)^(1,1) 


(7,12) 



continued on next page . . . 
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Table S4a 



continued from previous page 



Reaction 


Catalyst 


(1,1) + (0,1) ^(1,2) 


(7,12) 


(1,2) + (1,2) ^(2, 4) 


(7,12) 


(1,1) + (2, 4) ^(3, 5) 


(7,12) 


(2, 4) + (3, 5) ^(5, 9) 


(7,12) 


(1,1) + (5, 9) ^(6, 10) 


(7,12) 


(1,2) + (6, 10) ^(7, 12) 


(7,12) 


Generation 6 of length 25 (k = 15000) 


(0,1) + (1,0)^(1,1) 


(12,13) 


(1,1) + (1,0)^(2,1) 


(12,13) 


(2,1) + (0,1) ^(2, 2) 


(12,13) 


(2, 2) + (0,1) ^(2, 3) 


(12,13) 


(1,0) + (2, 3) ^(3, 3) 


(12,13) 


(2,1) + (3, 3) ^(5, 4) 


(12,13) 


(2, 3) + (5, 4) ^(7, 7) 


(12,13) 


(7, 7) + (2, 2) ^(9, 9) 


(12,13) 


(7, 7) + (3, 3) ^(10, 10) 


(12,13) 


(9,9) + (2,3)^ (11,12) 


(12,13) 


(11, 12) + (1,1) ^(12, 13) 


(12,13) 


Generation 7 of length 30 (k = 20000) 


(0,1) + (1,0)^(1,1) 


(13,17) 


(1,1) + (0,1) ^(1,2) 


(13,17) 


(1,2) + (1,1) ^(2, 3) 


(13,17) 


(1,0) + (2, 3) ^(3, 3) 


(13,17) 


(3, 3) + (1,1) ^(4, 4) 


(13,17) 


(3, 3) + (1,2) ^(4, 5) 


(13,17) 


(2, 3) + (4, 5) ^(6, 8) 


(13,17) 


(2, 3) + (6, 8) ^(8, 11) 


(13,17) 


(6, 8) + (4, 4) ^(10, 12) 


(13,17) 


(10, 12) + (2, 3) ^ (12,15) 


(13,17) 


(12, 15) + (1,2) ^(13, 17) 


(13,17) 


Generation 8 of length 35 (k = 27000) 


(1,0) + (0,1)^(1,1) 


(15,20) 


(0,1) + (1,1) ^(1,2) 


(15,20) 


(1,1) + (1,2) ^(2, 3) 


(15,20) 


(1,2) + (2, 3) ^(3, 5) 


(15,20) 


(3, 5) + (1,1) ^(4, 6) 


(15,20) 


(1,1) + (4, 6) ^(5, 7) 


(15,20) 


(3, 5) + (4, 6) ^(7, 11) 


(15,20) 



continued on next page . . . 
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Table S4a 



continued from previous page 



Reaction 


Catalyst 


(7, 11) + (4, 6) ^(11, 17) 


(15,20) 


(1,0) + (11, 17) ^(12, 17) 


(15,20) 


(1,1) + (12, 17) ^(13, 18) 


(15,20) 


(13, 18) + (1,2) ^(14, 20) 


(15,20) 


(1,0) + (14, 20) ^(15, 20) 


(15,20) 


Generation 9 of length 40 (k = 35000) 


(0,1) + (0,1) ^(0,2) 


(14, 26) 


(1,0) + (0,2) ^(1,2) 


(14, 26) 


(1,2) + (1,2) ^(2,4) 


(14, 26) 


(1,2) + (2, 4) ^(3, 6) 


(14, 26) 


(3, 6) + (1,0) ^(4, 6) 


(14, 26) 


(3, 6) + (2, 4) ^(5, 10) 


(14, 26) 


(5, 10) + (1,2) ^(6, 12) 


(14, 26) 


(6, 12) + (2, 4) ^(8, 16) 


(14, 26) 


(2, 4) + (8, 16);=^ (10, 20) 


(14, 26) 


(0,1) + (10, 20) ^(10, 21) 


(14,26) 


(10,21) + (0,2) ^ (10,23) 


(14, 26) 


(0,2) + (10, 23) ^(10, 25) 


(14,26) 


(10, 23) + (1,2) ^(11, 25) 


(14, 26) 


(8, 16) + (5, 10) ^ (13,26) 


(14, 26) 


(4, 6) + (10, 20) ^(14, 26) 


(14,26) 


Generation 10 of length 45 (k = 50000) 


(0,1) + (1,0);:^ (1,1) 


(18,27) 


(0,1) + (1,1) ^(1,2) 


(18,27) 


(1,2) + (1,0);=^ (2, 2) 


(18,27) 


(1,2) + (1,1) ^(2, 3) 


(18,27) 


(1,2) + (2, 2) ^(3, 4) 


(18,27) 


(2, 3) + (2, 3) ^(4, 6) 


(18,27) 


(4, 6) + (1,1) ^(5, 7) 


(18,27) 


(4, 6) + (3, 4) ^(7, 10) 


(18,27) 


(7,10) + (7,10) ^ (14,20) 


(18,27) 


(0,1) + (14, 20) ^(14, 21) 


(18,27) 


(14, 21) + (2, 2);^ (16, 23) 


(18,27) 


(2, 3) + (16, 23) ^(18, 26) 


(18,27) 


(14,21) + (4,6) ^ (18,27) 


(18,27) 
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Table S4b: List of reactions in the catalyzed chemistry and their catalysts. The table lists all 
the reactions that are part of the catalyzed chemistry with all its catalysts. The catalysts that belong 
to different generations (Gi to Giq) have been separated in different columns. It is easy to see from this 
table the amount of overlap between any two nested ACSs. For example, between the ACSs G5 and Gq 
which contain, respectively, 7 and 11 reactions, only one is common. 



Reaction 


Catalyst 


Gi 


G2 




G4 


G5 


Ge 


Gr 


Gs 


G9 


Gio 


(0,1) + (0,1) (0,2) 








(7,8) 










(14,26) 




(1,0) + (0,1) (1,1) 






(4,6) 




(7,12) 


(12,13) 


(13,17) 


(15,20) 




(18,27) 


(1,0) + (1,0) (2,0) 


(3,0) 


(4,2) 


















(0,2) + (1,0) (1,2) 








(7,8) 










(14,26) 




(1,1) + (0, 1) (1,2) 










(7,12) 




(13,17) 


(15,20) 




(18,27) 


(0,1) + (2,0) (2,1) 




(4,2) 


















(1, 1) + (1,0) (2, 1) 






(4,6) 






(12,13) 










(2,0) + (1,0) (3,0) 


(3,0) 




















(2,1) + (0, 1) (2,2) 




(4,2) 








(12,13) 










(1,1) + (1, 1) (2,2) 






(4,6) 
















(1,2) + (1,0) (2,2) 








(7,8) 












(18,27) 


(0, 1) + (2,2) (2,3) 






(4,6) 






(12,13) 










(1,2) + (1, 1) (2,3) 














(13,17) 


(15,20) 




(18,27) 


(1,2) + (1,2) (2,4) 










(7,12) 








(14,26) 




(1,0) + (2,3) (3,3) 












(12,13) 


(13,17) 








(2, 1) + (2, 1) (4,2) 




(4,2) 


















(1,2) + (2,2) (3,4) 




















(18,27) 


(1,1) + (2,4) (3,5) 










(7,12) 












(1,2) + (2,3) (3,5) 
















(15,20) 






(2,2) + (2,2) (4,4) 








(7,8) 














(3,3) + (1,1) (4,4) 














(13,17) 








(1,2) + (2,4) (3,6) 


















(14,26) 




(3,3) + (1,2) (4,5) 














(13,17) 








(2,1) + (3,3) (5,4) 












(12,13) 










(2,3) + (2,3) (4,6) 






(4,6) 














(18,27) 


(3,5) + (1,1) (4,6) 
















(15,20) 






(3,6) + (1,0) (4,6) 


















(14,26) 




(1,2) + (4,4) (5,6) 








(7,8) 














(1,1) + (4,6) (5,7) 
















(15,20) 




(18,27) 


(2,4) + (3,5) (5,9) 










(7,12) 












(5,6) + (1,2) (6,8) 








(7,8) 














(2,3) + (4,5) (6,8) 














(13,17) 








(2,3) + (5,4) (7, 7) 












(12,13) 










(3, 6) + (2, 4) ^(5, 10) 


















(14,26) 




(5, 6) + (2, 2) ^(7, 8) 








(7,8) 














(1,1) + (5, 9) ^(6, 10) 










(7,12) 












(4, 6) + (3, 4) ^(7, 10) 




















(18,27) 



continued on next page . . . 
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Table S4b 



continued from previous page 



Reaction 


Catalyst 


Gi 


G2 




G4 


Gs 


Ge 


Gt 


Gg 


Gg 


Gio 


(5, 10) + (1,2)^^(6,12) 


















(14,26) 




(3, 5) + (4, 6) -(7, 11) 
















(15,20) 






(7, 7) + (2, 2) (9, 9) 












(12,13) 










(1,2) + (6, 10) (7, 12) 










(7,12) 












(2, 3) + (6, 8)^^(8,11) 














(13,17) 








(7,7) + (3,3) ^ (10, 10) 












(12,13) 










(6,8) + (4,4) ^ (10, 12) 














(13,17) 








(9,9) + (2,3) - (11,12) 












(12,13) 










(6, 12) + (2, 4) ^(8, 16) 


















(14,26) 




(11,12) + (1,1) ^ (12,13) 












(12,13) 










(10,12) + (2,3) - (12,15) 














(13,17) 








(7, 11) + (4, 6)^^(11,17) 
















(15,20) 






(1,0) + (11, 17) -(12,17) 
















(15,20) 






(2, 4) + (8, 16) - (10,20) 


















(14,26) 




(12,15) + (1,2) ^ (13,17) 














(13,17) 








(0,1) + (10, 20) (10,21) 


















(14,26) 




(1,1) + (12, 17) ^(13,18) 
















(15,20) 






(10,21) + (0,2) ^ (10,23) 


















(14,26) 




(13, 18) + (1,2) ^(14,20) 
















(15,20) 






(7, 10) + (7, 10) (14,20) 




















(18,27) 


(0,2) + (10,23) (10,25) 


















(14,26) 




(0,1) + (14, 20) ^(14,21) 




















(18,27) 


(1,0) + (14, 20) ^(15,20) 
















(15,20) 






(10, 23) + (1,2) ^ (11,25) 


















(14,26) 




(8, 16) + (5, 10)^^(13, 26) 


















(14,26) 




(14,21) + (2,2) (16,23) 




















(18,27) 


(4, 6) + (10, 20) ^(14,26) 


















(14,26) 




(2, 3) + (16, 23) (18,26) 




















(18,27) 


(14, 21) + (4, 6) ^ (18,27) 




















(18,27) 
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Figure 1: Concentrations in uncatalyzed chemistries with a single food source. (A) Evolution 
of concentrations with time for a chemistry with kf = kr ~ A ~ (f> = 1. For simulation purposes, the size 
of the largest molecule was taken to be iV = 100. (B) Steady state concentration as a function of molecule 
size. Parameters take the same values as in (A) except that four values of (j> are shown, = 0, 0.1, 1, 10. 
Inset shows the same on a semi-log plot; the straight lines are evidence of exponential damping of Xn for 
large n (Eq. ([31)), with A = 1,0.77,0.49,0.16 for the four cases, respectively A is computed from the 
slope of a straight line fit after ignoring the smaller molecules (up to n = 4 in this case). 




Figure 2: Steady state concentration profile in an uncatalyzed chemistry with / = 2. The 

3D plot shows the concentration Xn of the molecule n = 712) as a function of ni and 712 in the steady 
state, for an uncatalyzed chemistry with kf = kr = a;(i^o) = a;(o.i) = 4> = I, N ~ 40. The inset shows a 
'top view' of the (711,712) plane with Xn indicated in a colour map on a logarithmic scale. 
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Figure 3: Steady state concentration profile for ACS65 (Eqs. (fS])). In all the cases kf = kr ~ 
A ~ 1,N ~ 100. (A) The concentration profile for two values of k for <f> = 15. (B) The concentration 
profile for four values of cj) for k = 3.0 x 10^. (C) The concentration profile for k = 2.5 x 10^, cf) = 15 
but with reaction (|5ap removed from the ACS (red curve) compared with the profile for the spontaneous 
chemistry, k = 0, (/> = 15 (green curve). The inset shows the same with x„ on a logarithmic scale. On 
the linear scale the two curves are indistinguishable. 




Figure 4: Steady state concentration profile for ACS(8,10) in a chemistry with / = 2. The 

molecules and reactions of the ACS are given in Supporting Table SI, the largest molecule being (8,10). 

kf ^ kr ^ X(i^o) = a;(o.i) = 1, = 10, k = 10^, N = 40. 
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Figure 5: Bistability in the dynamics of ACS65. 'Hysteresis curve' of the steady state concentration 
of A(65) versus k for kj ~ kr ~ A ^ 1,4> ^ 15,N ^ 100. The curve is obtained by using two different 
initial conditions (i) the standard initial condition x„ ~ for all n > 2, and (ii) a 'high' initial condition 
= 1 for all n > 2. In region 1 {k < = 6617) both initial conditions lead to a single fixed point in 
which X65 is very low, 10~^°. In region III {k > k^^ = 2226342) both initial conditions again lead to a 
single fixed point but in this fixed point xq^ is high, close to unity. In region II (k^ < k < k^^) the initial 
condition (i) leads to the lower fixed point and (ii) leads to the upper one. The transitions are very sharp, 
e.g., at K = 2226341 the system is numerically clearly seen in region II and at 2226343 in region III. 
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Figure 6: Phase diagram and concentration profiles for ACS4. (A) The steady state concen- 
tration versus k for kf = kr = A ~ 1, (/) = 25, = 15. The bistable region exists for the range 
< K < in which different initial conditions lead to two distinct steady state values of X4. The solid 
black curves correspond to the two stable fixed points, and the dotted black curve to the unstable fixed 
point. (B) The dependence of (red curve) and k^^ (blue curve) on cj) ioi kf = kr = A = 1, N = 15. 
The bistable region lies between the two curves; in the rest of the phase space the system has a single 
fixed point. The inset shows the location of the critical point {R, (j)); there is no bistability for (p < 4>. (C) 
Dependence of the phase boundaries on kf ioi kj. = A = 1, N = 15, with the inset showing the behaviour 
on a log-log plot. (D) The steady state concentration profile of molecules shown at five representative 
points in the phase space (numbered 1 through 5 and marked in (B)). Note that at the phase point 2 that 
lies between the and k^^ curves there are two steady state profiles corresponding to the two stable 
fixed points of the system. The figure marked 2a shows the profile starting from the standard initial 
condition, and 2b from the initial condition where a;„ = 1 for all n. 
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Figure 7: The dependence of the bistable region on catalyst length L. (A) The dependence 
of on L. (B) The dependence of on L. Simulations were done for extremal ACSs of length L 
generated by three algorithms (see Methods), represented in the figure by different colours. For each L 
the ACS in question has the property that the largest molecule produced in the ACS has L monomers 
and catalyzes all the reactions in the ACS. All simulations were done for kf — kr~A~(f>~l. N ~ 100 
in all cases except the curve for Algorithm 1, where N = 200, because in this case 'finite- iV' effects 
were quite significant at N ~ 100. The figures suggest an approximate power law growth of and 
exponential growth of k^^ with L. 




Figure 8: Reinforcement of a larger ACS by a smaller one: The case of ACS3+8. The figure 
shows the steady state concentration xs (in colour coding as indicated) for two different initial conditions 
as a function of K3 and ks, the catalytic strengths of A(3) and A(8) respectively. All simulations were 
done ioT k f ~ kr ^ A = 1, (j) = 20, N = 100. The two figures (A) and (B) differ in the initial condition of 
the dynamics. (A) The standard initial condition, (B) initial condition x„ = 1 for all n = 2, 3, . . . , A'^. 
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Figure 9: The product k xl is of order unity. This figure is produced from tlie same data as was 
used for Fig. [T] Simulations were done for chemistries containing extremal ACSs of length L generated 
by the three algorithms discussed earlier, represented in the figure by different colours. For this figure 
each chemistry was simulated at a value of k equal to the initiation threshold k^^ corresponding to that 
chemistry, and the steady state concentration xl oi the catalyst was determined in the low fixed point 
(starting from the standard initial condition). The parameters values are the same as in Fig. [T] (A) The 
product of and xl as a function of L. (B) xl versus k^^ on a log-log plot. The slopes of the fitted 
straight lines vary in the range -1.13 to -1.16 for the three algorithms (slope = -1 would have meant that 
K^^X2 is strictly constant. The figure shows that while each individual factor k^^ and x^ ranges over 
several orders of magnitude, their product, though not constant, is of order unity. 

































































IE 








- c 


■IE 


C 

1 


] - 































B 



CM 

c 











































































LJ 

ra. : 






■H-EilE 
" 1 


1 



















c 























































































- 1 


...r 












rl-i 


< 









01 2345678 



01 2345678 



01 2345678 



Figure 10: Examples of nested ACS pairs with different degrees of overlap for / = 2. hi the 

three cases the reaction sets have (A) maximal overlap, (B) partial overlap, (C) no overlap. The blue 
and red squares marking the grid points indicate the identity of molecules produced in the two ACSs; 
blue filled squares correspond to the products of the smaller ACS, red unfilled squares to those of the 
larger ACS. The x and y axes denote the number of monomers of type (1,0) and (0,1), respectively, in 
the molecules. The green rhombuses represent the two monomers. 
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Figure 11: The dominance of a cascade of nested ACSs with a molecule of size 441 (ACS441). 

The molecules and reactions of this ACS are listed in Supporting Table S2. The red curves show the 
steady state concentration x„ of all the molecules as a function of their size n, starting from the standard 
initial condition; blue dots show the concentrations of the ACS molecules. Insets show the same on a 
semi-log plot. It is evident that the large ACS molecules acquire a significant concentration. The catalytic 
strengths of the ACS molecules depend upon their size n according to K(n) = 500 x n^-^, and for all cases 
kf = kr = 1, (j) = 50, N = 800. The three figures differ in the level of sparscncss of the spontaneous 
chemistry in which the ACS is embedded. The spontaneous chemistry in (A) is fully connected, in (B) 
has degree 20, and in (C) has degree 2. 
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Figure 12: Dominance of a cascade of nested ACSs with length 64 (ACS(36,28)) in a / = 2 
chemistry. (A) 3D plot showing the steady state concentration Xn of the molecule n = (ni,n2) as a 
function of ni and n2, starting from the standard initial condition. The colour coding is on a logarithmic 
scale of the concentration. (B) A 'top view' of the same so that the ACS molecules and background are 
more clearly distinguished. The colour coding here is on a linear scale of concentration. The food set and 
ACS molecules have the highest concentrations and stand out as black dots. The catalytic strengths of 
the ACS molecules depend upon their size L = ni+n2 according to k{L) = 500 x L^-^, and fcj = fc^ = 1, 
(f> = 10, N = 65. The spontaneous chemistry has degree 20. 
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Figure 13: Dominance of a cascade of partially overlapping nested ACSs (ACS(18,27)). (A) 

Steady state concentration profile starting from the standard initial condition. (B) Top view of the same. 
(C) Sequence of steady state concentration profiles as each successive ACS is added to the chemistry. 
The legend for (A) is the same as for Fig. [12K and for (B) and (C) the same as Fig. [T2j3. kf = /s,. = 1, 
= 15, A = 50, and the spontaneous chemistry has degree 5. 
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Figure S2.1: Dependence of A on N for uncatalyzed chemistries. A = kr = I. We determine A 
at different values of N by fitting the steady state profiles to an exponential, Xn = ce^'^" (excluding the 
concentrations xi to X4). We see a dependence of A on iV for values of < 50, but A becomes essentially 
independent of N for > 100. 
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Figure S2.2: Dependence of steady state concentrations on A^ for a chemistry that includes 

an ACS. The figure shows the steady state concentrations for the chemistry that includes ACS65 (Eq. 
5 in main text) for A^ = 100, 200, 300, 500. In all cases A = kf = kr ~ l,(f) = 5. The steady state concen- 
trations show an A^-dcpendence upto A^ 100, but for A^ > 200 the profile becomes A^-independent. 
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Figure S3.1: Behaviour of A as a function of A:/ and 0. The figure shows the dependence of A on 
kf and (j) for an uncatalyzed chemistry, keeping A ~ kr = 1, N = 100. The curved surface was made with 
parameter values in the range 0<A:/<2,0.01<(/)<2. A is found to be a monotonically increasing 
function of kf and a monotonically decreasing function of cj). For (j) — 0, there is an analytical solution 
A ~ kf (see Eq. (4) in the main text). This was verified numerically in the range < kf < I (see solid 
line at (/) = 0). In the region kf > 1 the numerical integration does not converge at = as the steady 
state solution (Eq. (4) of the main text) x„ ~ AA"^^ = "^^^f~^ numerically very large for large n. 
The dotted extension of the line (1 < fc/ < 2) is simply the analytical result. Note that for most of the 
phase-space A < 1, except for very small values of (f>. 




Figure S3. 2: Behaviour of A as a function of kf and kr. The figure shows the dependence of A on 
kf and fc^ for an uncatalyzed chemistry, keeping A — (f) — 1, N = 100. A is found to be a monotonically 
increasing function of kf and a monotonically decreasing function of kj. . 



49 



